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ABSTRACT 

The recent star formation (SF) in the early-type spiral galaxy M81 is characterized using imaging 
observations from the far- ultraviolet (UV) to the far- infrared (IR). We compare these data with models 
of the stellar, gas, and dust emission for sub-galactic regions. Our results suggest the existence of 
a diffuse dust emission not directly linked to the recent SF. We find a radial decrease of the dust 
temperature and dust mass density, and in the attenuation of the stellar light. The IR emission 
in M81 can be modeled with three components: 1) cold dust with a temperature < T c >= 18 ± 
2 K, concentrated near the HII regions but also presenting a diffuse distribution; 2) warm dust 
with < T w >= 53 ± 7 K, directly linked with the HII regions; and 3) aromatic molecules, with 
diffuse morphology peaking around the HII regions. We derive several relationships to obtain total IR 
luminosities from IR monochromatic fluxes, and we compare five different star formation rate (SFR) 
estimators for HII regions in M81 and M51: the UV, Ho, and three estimators based on Spitzer 
data. We find that the Ha luminosity absorbed by dust correlates tightly with the 24 /im emission. 
The correlation with the total IR luminosity is not as good. Important variations from galaxy to 
galaxy are found when estimating the total SFR with the 24 /jm or the total IR emission alone. The 
most reliable estimations of the total SFRs are obtained by combining the Ha emission (or the UV) 
and an IR luminosity (especially the 24 fim emission), which probe the unobscured and obscured 
SF, respectively. For the entire M81 galaxy, about 50% of the total SF is obscured by dust. The 
percentage of obscured SF ranges from 60% in the inner regions of the galaxy to 30% in the outer 
zones. 

Subject headings: dust, extinction — galaxies: individual (M81, NGC3031) — galaxies: photometry — 
galaxies: ISM — galaxies: spiral — galaxies: stellar content — infrared: galaxies 



1. INTRODUCTION 

The multiwavelength analysis of the stellar populations 
in nearby galaxies is one of the most powerful ways of 
understanding how galaxies form and evolve (see, e.g., 
[ Meurer et alJll99ft IConselke et al.ll2000t Ic^alzett|l2001t 
lP£rez^ConzalezetHdTl2003a^^ . 
Indeed, relatively bright and close objects are the best 
laboratories to study galaxy evolution, given that we can 
obtain high quality data over a wide range of wavelengths 
with incomparable depth and physical resolution .The 
detailed characterization of nearby galaxies is essential 
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to establish the benchmark to which analysis of galaxies 
at intermediate and high redshift must refer. 

The multiwavelength data can be used to relate the 
scale where star formation normally occurs (i.e., the 
typical size of an HII reg ion, 20-200 pc; see, e.g., 
lYouneblood fc Hunterlll993) to the scales of the dynam- 
ical structures (bulge, arms, bars) and the entire galaxy 
(i.e., several kpc). Moreover, the combination of data 
from the ultraviolet (UV) to the far-infrared (FIR) and 
radio also allows simultaneous analysis of the three main 
(emitting) components of a galaxy: the stars, the dust, 
and the gas. 

One of the main problems in the study of the stellar 
populations in unresolved galaxies is the determination 
of the dust attenuation of the stellar emission. Indeed, 
the correction for dust extinction is the main source of 
uncertainty in the estimation of important characteristics 
of entire galaxies (i.e., of integrated stellar populations), 
such as the star formation rate (SFR), age or m etall icity 
<IKenmcutUll99cl ICordon et alMM ICalzetlJl27)0l iBell 
2002). This problem is even more severe when study- 
ing galaxies with current star formation, given that this 
process occurs in zones where the gas and dust densities 
are high (and the higher the densities, the more intense 
the star formation), resulting in high extinctions with 
large spatial variations. This effect translates into signif- 
icant uncertainties concerning t he evolution of galaxies 
throughout the Hubble time fsee lHonkinsll2004l and ref- 
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erences therein). 

One way to overcome the problem of dust attenuation 
is to combine data covering the widest wavelength range 
possible. This range must include the UV, where the 
emission comes mainly from hot young massive stars, 
and the optical and near-infrared, which are normally 
dominated by more evolved low mass stars. Spectro- 
scopic observations and narrow-band imaging concen- 
trating on emission-lines are also useful to study the 
stellar populations (for example, the Ha emission is di- 
rectly linked to the youngest stars). Much work char- 
acterizing the stellar populations in galaxies using a va- 
riety of such datasets and technical app roaches has ap- 
peared in the literature (among others. iBell fc de Jond 
| 200lUPapovich et al.ll200ll IPerez-Gonzalez ; et aTl F2003c: 
Kauffmann et al.ll2003t IHeavens et al.ll2004jh 

A complementary way to study the star formation 
(especially the most recent) in galaxies is using mid- 
and far-infrared observations. The emission at these 
wavelengths c omes mainly from dust heated by stellar 
light (see. e.g..|W alterbo s fc Greenawahll996llKennicuttl 
119981 iDevriendt et alJll999|) . The total IR emission (in- 
tegrated from 8 to 1000 /im) has been claimed to be an 
optimum SFR tracer for entire galaxies (at least for the 
ones dominated by star formation and with high dust 
content), g iven that it is not affected by dust atten- 
uation (see ISauvage fc Thuanl 119921 ISanders fe Mirabell 
119961 iKennicuttl 119981 and references therein). The 
IR has been widely used for stu dying local and dis- 
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tant star-form i ng g a laxies (e.g., 
Saunders et alJ Il990t iFlores et abl 
1999; |Fraji£e^chinLet aD 12001J; Rowan- Robinson! 
Chary fc Elbazl20011IKewlev et al.l2002HGoldader et al 
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20021 IDale fc Heloul 120021: IPerez-Gonzalez et all 12005). 
Although it is generally assumed that the heating of the 
dust emitting in the IR is dominated by light coming 
from the young hot stars (which would mean that the 
IR is a robust SFR estimator), it remains unclear how 
much the older stellar populations contribute (i.e., how 
we should correct the integrated IR emission to obtain 
robust values of the current SFR). Moreover, the validity 
of the equations correlating the SFR and the IR emission, 
which were derived for entire galaxies, is not well estab- 
lished when working at sub-galactic scales (for example, 
with HII regions). Finally, another interesting topic un- 
der study is the role of the aromatic molecules (loosely 
termed polycyclic aromatic hyd rocarbons, PAHs) in 
the star formation event (see iLeger fc Pugeij 11984 
lAllamandola et all Il98l iPuget fc Leger 119891 and 



most recently. [F^^j^& hreiber etHuT2004 : lBoseni et al.l 
l200llPeeters et al.H200|) ~ 

The study of the stellar populations in galaxies, and 
more specifically of the most recent star formation, has 
experienced remarkable progress in the past few years, 
as we have been able to gather observations from the 
UV to the FIR with similar depths and spatial resolu- 
tions. Two space missions have contributed significantly 
to this observation al advance: the G alaxy Evolution 
Explorer (GALE X: iMartin et all 120051) and the Spitzer 
Space Telescope (jWerner et al.ll2004a|K which have un- 
precedented sensitivity and angular resolution and cov- 
erage in the UV and IR, respectively. These facilities, 
jointly with other ground-based and space telescopes, are 
being used by the Spitzer Infrared Nearby Galaxy Sur- 



vey (SINGS, iKennicutt et al.ll2003j) to collect very high 
quality data from the UV to the FIR for a large sam- 
ple of nearby resolved galaxies representative of the local 
galaxy population. The main goal of SINGS is to provide 
a complete multiwavelength dataset to understand the 
formation of stars in galaxies and the properties of the 
dust. The SINGS dataset, including up to 7 broad-band 
measurements in the IR, and extensive ancillary data at 
optical and radio wavelengths, has prov ed to be very use- 
ful for exploring the se topics (see, e.g.. IDale et al.ll2005l 
iCalzetti et alll2005|) . 

In this paper, our main goal is to study the character- 
istics of the most recent star formation in a galaxy, test- 
ing the IR emission as a SFR estimator for sub-galactic 
regions by empirically analyzing its correlation with the 
"classical" star formation tracers, the UV continuum and 
the ionized hydrogen emission-lines (in our case, the Ha 
emissio n) . This paper c omple ments a similar analysis of 
M51 by Calzetti et al. (2005). Combining the data for 
the three SFR estimators, we can account for: 1) the 
photons coming from the newly-formed stars which do 
not interact with anything or are just scattered, and we 
detect in the UV; 2) the photons that interact with the 
gas, and are reemitted in emission lines; and 3) the pho- 
tons that interact with the dust, and are reemitted in the 
IR. By checking the consistency among these indicators, 
a clearer picture of the star formation event is obtained. 
The results in this paper will be used to study the differ- 
ent stellar populations (i.e., the star formation history) 
of M81 in a future work (Perez-Gonzalez et al. 2006, in 
preparation). 

This work concentrates on one of the galaxies in the 
SINGS sample, M81 (NGC3031), the largest in angu- 
lar extent in the survey (with a si ze of 27' x 14'). M81 
is a f ace-on spiral [SA(s)ab type, Ide Vaucouleurs et al 
1992| at a distance of 3.63±0.3 4 Mpc l|Freedman et al 
1994; iKarachentsev et alJ 12002') . At this distance, the 
resolution of the MIPS images corresponds to 100- 
700 pc (depending on the channel). This resolution 
allows the study of the star formation in M81 at dif- 
ferent scales ranging from typical sizes of HII regions 
to the size of large star-forming complexes and dy- 
namical structures (bulge, arms). M81 is an object 
with moderate global star formation (SFR = 0.4 — 
0.8 A^yr" 1 : iDevereux et alJ Il995l ILin et all l200l 
Gordon et alJ 12004) and extinction (MV) = 0.5 - 1.0, 
Garnett fc Shields! 119871 iKong et al J 12000) . It forms 



part of a group of about 25 galaxies with clear interac- 
tions among t hem (for example, am ong M81, M82, and 
NGC3077; see ISolinger'eT^IT977Tl. M81 also harbors 
an active galactic nucleus llPeimbert fc Torres-Peimbertl 
119811 iFilippenko fc Sargentl|1988D classified as a L INER 
Svl.5-L8|fHeckmanl 1198(1 IShuder fc Osterbrockl fl98ll 
IHoetal.1 119971). The nucleus i s a strong x- 
ray ( [ E lvis _fe van Spevbroedd Il98l iFabbianol J198& 
iPetre et al.ll993HPage et all2003ULa Parola et all2004l) 
and U V i|Wu et al.lll980l iHill et al.lll992t iReichen et alJ 
1994) emitter. M81 has been extensively studied at all 
wavelengths from x-rays to the radio. In this work, we 
will make use of the UV and emission-line data pub- 
lished for some of the individual HII regions in M81 
(lOarnett fc ShieldsH987llKaufman et al.ll987tlHill et alJ 
119951 ILin et al.ll2003(l . jointly with the new Spitzer and 
GALEX data. 
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This paper is organized as follows. Section [5] describes 
the data compiled for this work. In Sectional we present 
the models of the stellar and dust thermal emission de- 
veloped to analyze the UV-to-FIR data of M81. In Sec- 
tion we discuss the contribution of stellar light to the 
mid-IR fluxes, and different ways to estimate the total 
IR luminosity from monochromatic fluxes measured by 
Spitzer. We also concentrate on the analysis of the dust 
properties in several sub-galactic regions in M8I. Sec- 
tion[S]compares different SFR estimators for selected HII 
regions. Section [S] summarizes our conclusions. 

2. THE DATA 

We have compiled a variety of imaging and spectro- 
scopic data for M81 covering the ultraviolet (UV), op- 
tical, near-IR (NIR), mid-IR (MIR), and far-IR (FIR) 
wavelengths. The sources for each piece of the dataset 
are described in the following subsections. 

2.1. The Spitzer data 

Observations of M81 in the MIR and FIR were car- 
ried out with Spitzer ( We rner e t al.| |20 04aT) as a par t 
of the SINGS Legacy Project ijKennicutt et aT] [2 0031. 
The Infrared Array Camera (1RAC. iFazio et all 12004(1 
obtained images of M81 at 3.6, 4.5, 5.8, and 8.0 /xm in 
May 2004. The data were processed with version S10.0.3 
of the Spitzer Science Center pipeline, and then mo- 
saicked to produce a single i mage for each chan nel with a 
scale of 0.76 arcsec pixel -1 (|Reean et alf2 004). The Full 
Width at Half Maximum (FWHM) of the Point Spread 
Function (PSF) is 1.7", 1.7", 1.9", and 2.0" for the 3.6, 
4.5, 5.8, and 8.0 /im channels, respectively (correspond- 
ing to a physical size of ^30-35 pc at the distance of 
M81). The IRAC images were calibrated using the fac- 
tors and aperture corrections (to properly measure the 
photo metry of extended sources) found in iReach et alJ 
(2005). These corrections range from 6% in channel 1 
to 26% in channel 4. The 1-a sensitivities of the images 
are [0.2, 2.2, 2.0, 1.9] arcsec" 2 at [3.6, 4.5, 5.8, 8.0] ^m 
([25.7, 23.0, 23.2, 23.2] mag(AB) arcsec" 2 ). 

Images at 24, 70 and 160 /im were taken with 
the Multiband Im aging Photometer for Spitzer (MIPS, 
iRieke et al.l 12004^1 at three different epochs (to remove 
artifacts) in November 2003 and October 2004. The 
three epochs were processed and mosaicked together us - 
ing the DAT reduction pipeline l)Oordon et alJ [20051. 
producing final mosaics with pixel scales half of the orig- 
inal values (1.25, 4.9, and 8 arcsec pixel" 1 for 24, 70, 
and 160 fim, respectively). These mosaics have a size 
of 1° x 2.5°, large enough to include the entire galaxy. 
The FWHM of the PSFs are 5.7", 16", 38" for the 24, 
70, and 160 /jm channels (corresponding to 100, 282, 
and 669 pc in M81). The 1-a sensitivities of the im- 
ages are [0.9, 7.0, 20.0] ,uJy arcsec" 2 at [24, 70, 160] /im 
([24.0,21.8,20.7] mag(AB) arcsec" 2 ). 

2.2. The GALEX data 

The Galaxy Evolution Explorer (GALEX) observed 
M81 on December 8 th 2003 for a total of 3089 seconds 
split across two orbits. Note that since GALEX only ob- 
serves during night-time the usable time per orbit is only 
~ 1700 seconds. Images in the two GALEX bands, far- 
ultraviolet (FUV; A e ff=151.6 nm) and near-ultraviolet 



(NUV; A e ff=226.7 nm), were obtained simultaneously. 
The UV light is separated in these two bands by a 
dichroic beam splitter that also acts as a field-aberration 
corrector. Individual UV photons are then detected us- 
ing micro-channel plates with crossed delay-line anodes. 
The photon lists created are processed by the GALEX 
pipeline to produce the corresponding intensity maps in 
cou nts per second with a fi nal scale of 1.5 arcsec pixel" 1 
fsee iMo rrissey et, alJ 12005(1 . The reader is referred to 
iMartin et all <|2005[1 for a more detailed description of 
the GALEX mission. 

The large field of view of the GALEX instrument (a 
circle of 1.2° in diameter) allowed it to cover the en- 
tire M81 galaxy, along with M82 and Holmberg IX, in 
one single pointing. The absolute flux calibration of 
the GALEX data is based on multiple observations of 
white dwarf standard stars. Zero-point errors are esti- 
mated to be of the order of 0.15 mag both in the FUV 
and NUV channels. The astrometry of the GALEX im- 
ages is based on the Tycho-2 catalog <|Hog et al.ll2000(l 
and it is found to have a RMS error less than 1 pixel. 
The FWHM of the PSF of the GALEX images is a 
function of the brightness of the source and the posi- 
tion on the detector. In the case of the observations of 
M81 the FWHM is approximately 5" and 6" (88 and 
106 pc at the M81 distance, respectively), respectively 
for the FUV and NUV bands. The 1-a sensitivities of 
the images of M81 arc [1.8, 1.7] x 10" 2 juJy arcsec" 2 
([28.3,28.3] mag(AB) arcsec" 2 ) for the FUV and NUV 
channels, respectively. 

2.3. Optical and near-infrared data 

We complement the UV GALEX data and the 
MIR/FIR Spitzer images with optical and NIR data ob- 
tained at different ground-based facilities. UBVRI ob- 
servations were carried out in 2004-2005 with the 2.3m 
Bok Telescope at Stew ard Observatory using the prime 
focus wide-field imager ( William s et al.ll2"004(l Standard 
reduction algorithms for CCD images were applied to the 
data. We took special care with the distortion correction 
of the large l°xl° images. The depths of these optical 
images are [7.1,7.5,11.0,12.0,2.6] x 10" 2 /xJy arcsec" 2 
([26.8,26.7,26.3,26.2,27.9] mag(AB) arcsec" 2 ) for the 
UBVRI bands (1-a sensitivity), respectively. 

NIR images (JHK S bands) were tak en from the 
Two Micron All Sky Survey (2MASS; Uarrett et all 
IIXXJ. The depths are [80,140,140] /zJy arcsec" 2 
([19.1,18.5,18.5] mag(AB) arcsec" 2 ) in JHK (1-a sen- 
sitivity), respectively. 

Finally, we have used th e narr ow-band Ha image 
taken from iGreenawalt et all l)1998H . The calibration of 
this image was carried out by measuring the Ha+[NII] 
fluxes for al l the HII regions cataloged in lKaufman et all 
(1987) and iLin et alJ ((2003') . Our calibration gives a 
total Ha+[NII] flux for M81 of 4.7 ± 1.0 ergs" 1 cm" 2 . 
T he difference between th is flux and the estimations 
of IGreenawalt et all (119981 4.8 ± 0.7 ergs" 1 cm" 2 ) and 
iDevereux et al.l l(1995l 4.53 ± 1.35 ergs" 1 cm" 2 ) is less 
than 4% and well within errors. In the following discus- 
sion, we will use pure Ha fluxes that have been measured 
from the narrow-band image and corrected for [Nil] con- 
tamination assuming an average [Nil] /Ha ratio of 0.4 
((Garnett fc Shields! Il987jl . The 1-a Ha sensitivity limit 
of this image is 1 x 10" erg s" 1 cm" 2 arcsec" 2 . 
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2.4. Photometry 

We produced two sets of images, one set composed 
of images convolved to match the PSF of the image in 
our dataset with the worst resolution, the MIPS 160 fim 
image, and another set convolved to the MIPS 24 /mi 
PSF, which is adequate to study the dust emission in 
individual HII regions. The World Coordinate System 
(WCS) for the UV, optical, NIR, and IRAC images was 
determined using 2MASS stars (more than 30 per frame) . 
For the MIPS images, given that stars are normally very 
faint or undetected, we determined the WCS by using 
bright point-like regions detected in the IRAC bands. 
The images were aligned and cropped to the same size 
using the WCS for each frame. 

For the image set convolved to the 160 /mi resolution 
(160RES set hereafter), we selected 49 circular regions in 
the arms corresponding to resolved peaks in the 160 /im 
emission. We also selected 6 regions in the zone between 
the arms and the nucleus dominated by a more diffuse 
emission, two regions around the nucleus (one enclosing 
the inner nucleus, and another the circumnuclear region), 
and an aperture enclosing the entire galaxy. The radii of 
the photometric apertures were chosen to enclose most of 
the emission from the peaks (ranging from 0.6 to 1.8 kpc). 
The photometric apertures are shown in the left panel of 
Figure n Most of the selected regions were identified in 
the 24 /im and Ha images with bright star-forming com- 
plexes, normally enclosing several individual HII regions 
(see right panel of Figure ^) . 

Proper background subtraction is essential when deal- 
ing with data spanning a wide range of wavelengths, es- 
pecially if some of the images include very bright diffuse 
em ission (such as the im ages of M81 in the MIR and FIR, 
see lGordon et alJl2004[) . We followed two approaches in 
the background estimation, each of them optimized for 
one of the two problems that we can address with this 
dataset. First, analyzing the properties of the global 
stellar population and the dust in several sub-galactic 
regions in M81 requires photometry that accounts for 
the entire emission of the region at any wavelength (in- 
cluding all stars, gas and dust in the region) , so the best 
approach is to estimate a global background, i.e., a back- 
ground that only includes the sky emission around the 
galaxy. However, second, for studying different SFR es- 
timators, one of the main goals of this paper, it is im- 
portant to account for the diffuse emission surrounding 
each star-forming complex. This emission is probably not 
directly related to the recent star formation, but arises 
from a source in the galaxy of a different origin (e.g., 
older stars, diffuse dust emission). 

As an example, the 160 /im image of M81 shows a 
bright diffuse component that is not linked with any star- 
forming region (detected in the Ha or 24 /mi image). 
Other authors have also detected an extended and diffuse 
dust component beyond the star-f orming disk in a vari- 
ety of galaxies (s e e, am o ng others. jHippelein et alJl2 00rt 
IPonescu fc Tuff3 l2003 lHm7et al J 120041 IStickel et a ,U 
120051 iPopescu et al.ll2005|) . The cold dust responsible 
for this emission is probably not (only) heated by the 
newly-formed stars, but by a more diffuse heating source, 
such as an evolved stellar population or Lya and Ly- 
man continuum pho tons escaping from HII regions (see 
IPopescu et afll200"l . 



The global background for M81 was measured in sev- 
eral zones located around the galaxy, which were ob- 
served at all wavelengths and clear of bright objects 
(green rectangles and elliptical annulus in the left panel 
of Figure The local background was measured in 
small circular apertures around the arms (green circles 
in the left panel of Figure ^) . For each individual re- 
gion studied, we averaged the values obtained for the 
two closest sky apertures, one in the inter-arm region 
and another in the outer part of the arm. Two photo- 
metric catalogs were produced using the two background 
estimations (marked as LOCALBKG and GLOBALBKG 
photometry) , and the implications of each choice are dis- 
cussed in the following Sections. 

For the image set convolved to the 24 /an resolution 
(24RES hereafter), we extrac ted photometry for 59 in- 
divid ual HII regions found in iGarnett fc Shields! (Q987) 
and iKaufman et "all l)1987h . These regions are shown in 
the right panel of Figure The original aperture sizes 
used by those authors were very small compared with 
the resolution of our 24 /mi image, so we extracted pho- 
tometry for circular apertures with a 7 pixels (~8.4") 
radius, which enclosed most of the emission for each HII 
region as seen in the Ha image (the radius of these aper- 
tures corresponds to 0.1 kpc), and is an adequate match 
to the 24 fim PSF. The sky estimation was performed 
locally in circular apertures near the individual HII re- 
gions 12 (green circles in the right panel of FigurcQJ. For 
the entire galaxy, the sky was measured in an elliptical 
annulus. 

Aperture corrections were applied to the fluxes of the 
regions selected in each image set to obtain the photome- 
try for an "infinite aperture" using theoretical STinyTim 
PSFs for the MIPS bands. All the fluxes were corrected 
for Galactic extincti on using a value of A(B) = 0.346 
(Schletffil et; alj|1998j) and the extinction law published 
bv lCardelli et alJ l|1989D with R v =3.1. The fluxes for 
all regions are given in Tables ^ an d 121 

3. MODELS OF THE UV-TO-FIR EMISSION OF M81 
3.1. The stellar and gas emission model 

For the regions selected in the 160RES image set, we 
built spectral energy distributions (SEDs) from 0.150 nm 
to 160 /mi. SEDs were also built for the HII regions 
selected in the 24RES image set from 0.150 /jm to 24 /im. 
We then fitted the SEDs with models of the stellar, gas, 
and dust emission. These models will be used in this 
paper to remove the stellar contribution to the emission 
in the MIR and FIR bands. In a future paper, the results 
on the stellar populations of M81 obtained with these 
models will be discussed in detail. 

The photometric points at wavelengths bluer than 
4.5 /mi were fitted to stellar population synthesis models. 
We also included the Ha equivalent width in the fits ( see 
iGil de Paz et alJl2000tlPerez-Gonzalez et al.ll2003bl(jh 

We assumed that the star formation history in each 
region can be described by a burst with an exponen- 
tially declining star formation rate with time scale r, 

12 Only one approach of the background determination was fol- 
lowed for the 24RES case (local background), given that we were 
interested in comparing the SFR estimators. This task benefits 
from a proper subtraction of the diffuse emission surrounding the 
HII regions. 
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TABLE 1 

Positions and photometry (GLOBALBKG and LOCALBKG cases) for the regions selected at 160 resolution. 



Field 


BKG 


M81 


Rc 


?02 




Reg04 


Rppf}5 




Reg07 






RA (J2000) 




09:55:32.2 


09:55:32.2 


09:55:32.2 


09:55:16.7 


09:55:05.0 


09:55:35.1 


09:55:53.3 


09:55:35.7 


09:55:26 


DEC (J2000) 




69:03:59.0 


69:03:59.0 


69:03:59.0 


69:06:12.4 


69:05:22.9 


69:06:24.1 


69:01:47.0 


69:07:29.6 


69:08:08 


radius 




680.0 


64.0 


104.0 


40.9 


50.8 


51.3 


56.0 


31.9 


31.6 


log [L(Ha)] 


G 


40.78 


39 


67 


39.40 


38.79 


38.88 


39.10 


38.86 


39.10 


39.24 


L 


40.69 


39 


53 


39.18 


38.40 


38.42 


38.82 


38.54 


39.06 


39.09 


log [£(3.6)] 


G 
L 


43.18 


42 


63 


42.38 


41.63 


41.62 


41.75 


41.81 


41.14 


41.10 


log [£(4.5)] 


G 
L 


42.84 


42 


28 


42.03 


41.27 


41.25 


41.40 


41.49 


40.81 


40.77 


log [£(5.8)] 


G 


42.63 


41 


98 


41.76 


41.05 


41.06 


41.21 


41.21 


40.73 


40.72 


L 


41.42 










39.30 


40.03 




40.32 


40.22 


log [£(8.0)] 


G 


42.61 


41 


76 


41.58 


40.93 


40.99 


41.18 


41.11 


40.88 


40.88 


L 


42.33 


40 


65 


40.79 


40.34 


40.42 


40.78 


40.58 


41.05 


40.70 


log [£(24)] 


G 


41.98 


41 


25 


40.90 


40.20 


40.29 


40.49 


40.38 


40.28 


40.35 


L 


41.92 


41 


13 


40.68 


39.90 


39.96 


40.31 


40.20 


40.27 


40.28 


log [£(70)] 


G 


42.67 


41 


83 


41.62 


40.98 


41.00 


41.23 


41.15 


40.98 


41.05 


L 


42.68 


41 


79 


41.52 


40.79 


40.72 


41.11 


41.08 


40.98 


41.00 


log [£(160)] 


G 


42.99 


41 


78 


41.78 


41.27 


41.37 


41.55 


41.46 


41.30 


41.31 


L 


42.99 


41 


67 


41.62 


41.05 


41.11 


41.38 


41.36 


41.29 


41.22 



Note. — Table 1 is available in its entirety via the link to the machine-readable version above. Units of right ascension are hours, 
minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. The radius of each region is given in arcsec. All 
luminosities refer to observed values (no internal or Galactic extinction corrections have been applied) and the units are erg s — 1 . The 
first row for each wavelength refers to the photometry measured with a GLOBALBKG estimation (second column set to the value G), 
and the second to the LOCALBKG case (second column set to the value L). 



TABLE 2 

Positions and photometry for the regions selected at 24 /im resolution. 





Field 


M81 


kauf240 


kaufl04 


kauflOl 


Munchl 


kaufl02 


kaufl98 


Nol7 


kaufl97 


RA (J2000) 


09:55:33.2 


09:56:27.0 


09:56:20.6 


09:56:19.9 


09:56:18.9 


09:56:18.3 


09:56:17.2 


09:56:17.2 


09:56:16.6 


DEC (J2000) 


69:03:55.1 


69:04:24.0 


69:01:12.1 


69:04:25.3 


68:49:40.9 


69:03:34.9 


69:06:07.4 


69:06:07.4 


69:05:36.2 


radius 


680.0 


8.4 


8.4 


8.4 


8.4 


8.4 


8.4 


8.4 


8.4 


A(V) 




1.05 


1.31 


1.01 


0.37 


1.91 


0.42 


1.60 


0.67 


log 


L(FUV)] 


42.74 


39.49 


39.67 


40.00 


39.51 


39.98 


40.22 


40.22 


39.88 


log 


L(NUV)] 


42.82 


39.25 


39.57 


39.93 


39.06 


39.83 


40.08 


40.08 


39.84 


log 


L(Ha)\ 


40.77 


37.79 


38.14 


38.24 


38.01 


38.41 


38.38 


38.38 


38.28 


log 


£(8.0)] 


42.61 


39.08 


39.64 


39.46 


38.50 


39.65 


39.64 


39.64 


39.64 


log 


£(24)] 


41.99 


38.73 


39.73 


39.00 


38.67 


39.35 


39.37 


39.37 


39.42 



Note. — Table 2 is available in its entirety via the link to the machine-readable version above. Units of right ascension are 
hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. The radius of each region is given 
in arcsec. All luminosities refer to observed values (no internal or Galactic extinction corrections have been applied) and the units 
are erg s . 



metallicity Z, age t, and attenuated by an amount 
described by the quantity A(V). The attenuation at 
any wavelength was calculated from t he free parameter 
A(V) using the lCharlot fc Falll pOOOj) recipe. The stel- 
lar e mission in our models was taken fr om the PEGASE 
code (|Fioc fc Rocca-Volmerangel I1997J) . and we added 
the emission from the gas heated by the stars (emission 
lines and nebular continuum) usi ng the emission and re- 
combination coefficients given bv lFerlanHI Ip8(l) for an 
ele ctron temperature T P — 10 4 K, the relations given 
bv iBrocklehurstJ l)1971|) . and the theoretical line-ratios 
expected for a low density gas (n e = 10 2 cm~ 3 ) with 
Tf. = 10 4 K in the recombination Case B l)Osterbrockl 
1989). We fixed the metallic ity of our models t o 
based on the results fro m iZar itskv et ail (119941) . 
iStauffer fc Bothunl (|1984l) . and lKong et alJ <|20n(l .^ho 
showed that the metallicity in M81 is nearly solar with 
a very weak radial dependence. We probed the solution 



space in the following ranges for the three remaining free 
parameters [t, t, A(V)]: 1) we assumed r values from an 
almost instantaneous burst (r — 1 Myr) to an almost 
constant SFR (r = 15 Gyr) using a logarithmic interval 
of O.ldex (in yr); 2) ages were probed from t = 1 Myr 
to t = 13 Gyr in logarithm intervals; 3) extinction val- 
ues ranged from A(V) — to A(V) = 3 in intervals of 
0.05 mag. 

3.2. The model of emission by dust 

After fitting the UV-to-NIR part of the SEDs with stel- 
lar/gas emission models, we subtracted the contribution 
from the stars/gas to the fluxes in the MIR and FIR. 
We then fitted models of dust emission to the remaining 
fluxes. These models are described by an equation of the 
form: 

F dust (A) = ]TcV(A)5 A (T,U) (1) 
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Fig. 1. — The regions for which we extracted photometry from M81 are shown superimposed on the MIPS 24 (im image. The left panel 
shows the effect of convolving to the 160 fim resolution, and the right panel show the photometric apertures for the HII regions selected 
at the original 24 /an resolution. On the left panel, blue circles enclose star-forming complexes in the arms. The red circles enclose regions 
of diffuse emission in the zone between the nucleus and the arms. The magenta region encloses the inner nucleus of M81. The area outside 
the magenta circle and inside the orange circles is defined as the circumnuclear region. Green rectangles show the zones used for estimating 
the global background (selected to be free of br ight sources in all the images). For the right panel, bl ue regions come from the sample of 
HII regions studied by Garnctt & Shields (1987), and red regions were found in Kaufma n et alj 119871) . For both panels, the black ellipse 
depicts the aperture used for the entire M81 galaxy. The green ellipse and circles show the zones where the local sky was measured (see 
text for details). 



where the sum extends over the number of dust com- 
ponents (each of them emitting as a modified black- 
body), C % = M^ust/D 2 , D is the luminosity distance 
to the source, M l dust are the masses of each dust com- 
ponent, and the wavelength dependent mass ab- 
sorption coefficients. We adopted a model with 3 compo- 
nents: 1) warm silicates; 2) cold silicates (both compo- 
nents with grains of sizes a ~ 0.1 /im); and 3) aromatic 
molecules. Mass absorption coefficients for astronomical 
silicates were co mputed from Mie theory using the dielec- 
tric functions of lLaor fc Drains! l)1993|) . C ross sections 
for th e aromatic molecules were taken from lLi fc Drairiel 
(120011). As the can onical aromatic spectrum (see, e.g., 
iWerner et al.ll2004b() exhibits no features beyond 20 /im, 
we re-computed the aromatic cross-sections (and mass 
absorption coefficie nts) leaving o ff the last three terms 
of Equation 11 in Li fc Draind (|2001^) . The aromatic 
component is included mainly for completeness. The 
aromatic emission is a stochastic process rather than 
an equilibrium process as assumed in Equation ^ This 
means that our model should underestimate the total 
aromatic masses. A proper treatment of the stochastic 
aromatic emission would require a full radiative transfer 
solution to derive the radiation field at each position in 
the galaxy. The attendant assumptions are beyond the 



scope of this paper, but will be addressed in forthcom- 
ing detailed models of individual regions. Our results are 
not affected by this issue, given that they only refer to 
the total dust masses, to which aromatic molecules make 
only a minor contribution. 

The fitting procedure is based on a Monte-Carlo tech- 
nique. With our assumptions about the components in 
the observed SED, we searched a limited set of param- 
eters for dust mass (through Ci) and temperature for 
each component. We assumed temperature ranges of 
200-750 K, 35-100 K, and 10-25 K for the aromatics, 
warm, and cold silicates, respectively. At each set of 
temperature points, a similar grid of mass scalings (Ci) 
was explored for each component. After an initial, coarse 
grid search, a finer temperature and mass scaling grid 
was defined for each component and the procedure re- 
peated. We iterated this procedure (typically only two 
refinements of the grid were needed) until a good fit was 
found based on a x 2 goodness of the fit estimator. While 
the results of this technique are not unique, it yields rep- 
resentative estimates of both the dust mass and temper- 
ature (and errors obtained through the analysis of the 
\ 2 distribution), given the input assumptions. Some ex- 
amples of the final fits to the entire SEDs are shown in 
Figure Note that our simplistic models (which only 
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0.1 1 10 100 1000 0.1 1 10 100 1000 



X [/J,m] X [jtim] 




X [/u.m] X [/U.m] 

Fig. 2. — Spectral energy distributions for the entire M81 galaxy (upper-left panel), the nucleus (upper-right), and two regions in the 
arms (lower panels) selected at 160RES, corresponding to the photometry estimated with the global sky subtraction. The open circles 
in the SED of the enti re galaxy refer to sev eral photometric points extracte d from NED (for apertures nearly as large as ours): three 
UV/optical fluxes from lCode fc WelcrJ J19821) andjde Vaucouleurs et alJ 119921) . the three NI R fluxes from 2MA SS IJarrett et al.120001). th e 
flux in the iV-band from Kleinmann fc Lowl 119701) . and the four IRAS fluxes measured by Rice et alJ {1988) and Devc reux et all 11993) . 
Horizontal error bars for each photometry point show the width of the used filter. The data have been fitted (red line) to a model of the 
emission of the stars (a single burst with a exponential SFR), the gas (including the most common emission lines, shown in the plot with 
a width equal to the resolution of the stellar model at each wavelength), and the dust (including two dust components of two different 
temperatures , and aromatic molecules). Our models of the dust emission are compared with the ones carried out following the method in 
Drainc et al. 12006, green line). 



assume two dust temperatures) usually produce an arti- 
ficial inflexion point at around 50-60 /im (where the cold 
and warm dust components join). This artifact is also 
present in other models (for example, in th e ones plotted 
in Figure El using the technique described in lDraine et al.l 
2006), but it is not seen in any IR spectrum. However, 
it does not affect our estimates of the integrated IR lu- 
minosities, given that the energy output is dominated by 
emission at longer wavelengths. 

The small number of independent observations in the 
MIR and FIR justifies the simplistic approach in our 
models with just three emitting components that need 
5 free parameters to fit. We will use these models to es- 
timate the integrated IR luminosities of regions in M81 
and to study the characteristics of the dust in those re- 



gions. We tested our results with other more complicated 
models (that use more free parameters and/or assump- 
tions) to check that our approach did not introduce any 
systematic effects in those areas. We use d comparison 
models developed bv iDale fcFIglgjj l (120021 see the next 
Section for more details) and lDraine et alJ l)2006l) to char- 
acterize all the SINGS galaxies. The latter models are 
based on the study of the Milky Wa y extinction carried 
out in iWe ineartncr & Drainc (2001). Briefly, these au- 
thors assume a certain stellar radiation field (based on 
the NIR data in our case) that heats the dust grains. The 
distribution of sizes of the dust grains is that observed 
for the Milky Way. The grain temperatures are calcu- 
lated using a complete treatment of temperature fluctu- 
ations following single-photon heating and by studying 
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the thermal equilibrium of the dust. Figure [21 shows the 
comparisons of our fits and the on es carried out with 
the models bv IDraine et all (|2006). The fits are very 
similar, showing large differences only for wavelengths 
longer than 160 /mi, where the models bv IDraine et alJ 
( 2006) sometimes overpredict the observed flux. Despite 
this, our estimations of the TIR luminosities are very 
similar to those obtained from the comparison models 
(the average difference is less than 4%, with a scatter of 
6%). However, the dust masses we obtain with our mod- 
els are systematically smaller (on average, 0.2±0.2 dex). 
This difference is linked to the sizes and composition of 
the grains assumed by each set of models. While we as- 
sume a fixed dust g rain size and two dust temperatures, 
IDraine et all 120061 assume a dust size distribution based 
on Milky Way abundances, producing a range of temper- 
atures. This size distribution includes larger amounts of 
cold dust than our two component fits. The difference 
in dust masses is comparable to the quoted uncertainties 
(see Figure El and does not affect our results. 

3.3. Degeneracies of the models 

Our stellar population synthesis models are subject to 
the typical degene racies in this ki nd of study, for exam- 
ple a ge-metallicitv (|Worthe v) 19941) or extinction-age (see, 
e.g.. |Gil de Paz fc Madorej |2002f iPerez- Gonzalez et al.l 
2003bl and lKauffmann et al.ll2003[) . 

The well-known age-metallicity degeneracy was arti- 
ficially removed in our models by fixing the metallic- 
ity to t he solar value based on estimates by several 
authors llZaritskv et all 11991: iStauffer fc Bothunl 119841: 
IKong et al.ll2000fl . 

The main source of degeneracy in our stellar emission 
models then comes from the lack of knowledge about 
the star formation history of the galaxy. Indeed, sim- 
ilar SEDs can be obtained with a young instantaneous 
(r ~ 0) burst or with older and less intense continu- 
ous star formation (r 3> 1). Therefore, our simplistic 
assumption of an exponentially decaying burst of star 
formation can be very far from the reality, where sev- 
eral small bursts may have occurred at different epochs. 
The final photometry would be, in general, blind to that 
formation history. However, given that M81 is an early- 
type spiral, it would be expected to have an old and 
massive burst forming the bulge, where most of the stars 
of the galaxy lie. This old stellar population would dom- 
inate the emission in the NIR. More recent star forma- 
tion (switching on and off during the life of the galaxy) 
would be concentrated in the disk, with the most recent 
events contributing significantly to (or even dominating) 
the UV emission. Extinction can make the two scenarios 
previously described nearly indistinguishable. 

However, the extinction and the lack of knowledge 
about the star formation history have only a minor effect 
on our results. In this paper, the models have been used 
to estimate the contribution of stellar light to the fluxes 
in the IRAC and MIPS bands. At these wavelengths, the 
extinction is negligible. Variations in the age or the star 
formation history do not strongly affect the shape of the 
SED in the NIR, particularly when the emission at these 
wavelengths is dominated by old stars (at least, for an 
early type spiral such as M81). 

In the case of the dust emission models, there are to 
main degeneracies. The first comes from the lack of 



knowledge about the aromatic spectrum in M81, and 
about the possibility of changes of that spectrum from 
region to region. Given this degeneracy, we will avoid 
any discussion about aromatic molecules in this paper. 
The second degeneracy is linked to the limited data in the 
FIR, which is unable to locate the emission peak accu- 
rately. Indeed, the SEDs seem to continue rising beyond 
our reddest point at 160 //m (or the maximum is reached 
at approximately that wavelen gth). Based o n our m od- 
els and the ones developed bv IDraine et al.l l)2006[) . we 
estimate that the lack of data at A > 160 /im leads to 
an uncertainty of less than a factor of 2 in the derived 
integrated IR luminosities and dust masses (which would 
be underestimated by our models if a significant amount 
of cold dust -with T < 15 K- is also present in M81). 
This uncertainty is smaller than the quoted errors. The 
absence of any data in the MIR between 24 and 70 /jm 
should not affect the results significantly, since the en- 
ergy in this wavelength regime is a small fraction of the 
total in the IR. 

4. PROPERTIES OF THE DUST IN M81 

In this Section, we will describe our main results on 
the properties of the dust in M81. All the results in this 
Section were obtained from the 160RES dataset because 
the dust emission is better characterized in this case (FIR 
data is available), and they span a larger range of radial 
distances, mapping most of the galaxy. We analyze the 
stellar and dust emission in the MIR from a statistical 
point of view, and estimate the integrated IR luminosity, 
which will be used in the comparison of SFR estimators 
in Section We will specifically discuss the results ob- 
tained with the LOCALBKG photometry, where the cal- 
culated integrated luminosities refer to the dust heated 
only by the young stars, since the more diffuse emission 
(of more uncertain origin) is removed. We also include a 
discussion about the GLOBALBKG photometry, which 
might be useful for the study of the global stellar content 
in sub-galactic regions. In the last part of this Section, 
we will concentrate on the radial dependence of the dust 
properties in M81. Given that we are interested in the 
global characteristics of the dust, the results in the last 
subsection refer to the photometry measured with the 
regions selected at 160RES where the background was 
estimated globally (to include all the emission). 

4.1. Stellar emission in the MIR 

Our stellar population models show that the emission 
of the regions selected at 160RES is completely domi- 
nated by stars at 3.6 /im and 4.5 fim, with very little 
contamination from the dust: the median and standard 
deviation of the fraction of the total emission coming 
from the dust for all the 160RES regions is ± 15% at 
3.6 /im and ± 10% at 4.5 /im. Stellar emission is also 
dominant at 5.8 /im (40 ± 20% of the flux comes from 
the dust). At redder wavelengths, most of the emission 
comes from the dust: 80 ± 20% of the emission at 8 /im, 
93 ±7% at 24 /im, and 100±1% at 70 and 160 /im. These 
fractions are not affected by the intrinsic degeneracies of 
the modeling procedure, given that the models in this 
part of the spectrum are not very sensitive to changes 
in age or metallicity (the SEDs are dominated by the 
Ray leigh- Jeans spectrum of cold old stars), and the NIR 
and MIR are not strongly affected by extinction. 
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The fractions of the dust emission in the IRAC bands 
were checked using the recipe found in iPahre et al.l 
(2004). The emission at 3.6 and 4.5 ^m is assumed to 
be completely dominated by stars and we can estimate 
the contribution of the stellar emission to the 5.8 and 
8.0 /im channels using the typical colors of an MO III 
star [m(3.6) — m(4.5) = —0.15 mag, m(4.5) — m(5.8) = 
+0.11 mag, m(5.8) - m(8.0) = +0.04 mag in the 
Vega system, corresponding to F l/ (3.6)/F u (i.5) = 1.78, 
^(4.5)7^(5.8) = 1.39, and F u {h.8)/F v (8.0) = 1.78], 
which should dominate the stellar emission at these wave- 
lengths. Using this recipe, we obtain an average dust 
contribution to the 5.8 and 8.0 /im emissions of 45 + 26% 
and 77 ± 22%, respectively. 

Note that the regions used at 160RES enclose both 
star- forming complexes and/or zones of the arms where 
the emission is most probably dominated by older stars. 
The size of the photometry apertures has an effect on the 
stellar emission fractions previously quoted: for actual 
HII regions (i.e., for smaller apertures such as the 24RES 
ones) where the emission from young stars and warm 
dust is relatively more important, the fractions change 
by significant factors. According to our stellar popula- 
tion models, the 3.6 /im and 4.5 /im emissions have a 
non-negligible dust component: 10 ± 10% and 10 ± 10% 
of the total 3.6 and 4.5 /im flux comes from the dust, re- 
spectively. Note that heavily extincted young stars that 
would be invisible in the optical could contribute to the 
observed flux in the NIR at 24RES. According to our 
models, the dust emission dominates the total flux at 
redder wavelengths: 70 ± 20% at 5.8 /im, 90 ± 10% at 
8.0 /urn, and 100 ± 10% of the total emission at 24 /im. 

4.2. Integrated IR emission 

In this Section, we use the regions selected in the 
160RES image set to investigate possible relationships 
between monochromatic fluxes and colors in the MIR 
and FIR and the total emission of the dust integrated 
from 8 to 1000 /im [TIR luminosity, L(8 - 1000) here- 
after]. The following discussion will be mainly based on 
the photometry obtained for a GLOBALBKG. We will 
analyze the effect on these results of the diffuse IR emis- 
sion by also using the photometry obtained with a LO- 
CALBKG. This should give relationships which are bet- 
ter suited for SFR studies. The contribution of the stellar 
light has been subtracted from the MIR-FIR photometry 
prior to determining L(8 — 1000). 

4.2.1. The estimation of L(8 - 1000) 

TIR luminosities were calculated for the 160RES re- 
gions by integrating the fitted dust models. We also 
calculated the integrated luminosity from 3 to 1100 /im 
[L (3 — 1100)] to compa re our data wit h the se mi-empirical 
models developed bv iDale fc Heloul l)2002f) . These au- 
thors obtained an equation (Equation 4 in that paper) 
to estimate the L(3 — 1100) from the three MIPS fluxes. 
This relationship is based on models of the global emis- 
sion of normal star-forming galaxies with L(3 — 1100) > 
10 8 Lq. Thus, the application to faint individual HII re- 
gions within a galaxy may not be justified. In fact, the 
SEDs for sub-galacti c regions in M8 1 are c older than the 
coldest model in the IDale &: Heloul l)2002|) template set, 
as also noted by IDale et aLl ~ ij2"005f) for sub-galactic re- 



gions in M51 and NGC7331, and bv lHelou et all lEffll 
for NGC300. 

For the GLOBALBKG photometry, our own estima- 
tions of L(3 — 1100) are systematically fainter (7% on 
average, with a scat ter of 6%) than the ones obtained 
using Equation 4 in IDale fc Heloul lj2002|) . The scatter 
is mainly due to the two nuclear regions and the coldest 
regions in the inter-arm zone; without these regions, the 
scatter is 2%. For the entire galaxy, our models give a 
value of L(3 — 1100) that i s 1% lower than the one ob- 
tained with the equation in IDale fc Heloul l)2002f) . These 
small systematic differences (comparable to the uncer- 
tainties in the measured MIPS fluxes) are consistent with 
the fact that the sub-galactic regions in M81 are colder 
than the coldest model in the IDale &: Heloul (|2002j) tem- 
plate set. 

If we use the photometry calculated with the LO- 
CALBKG, the values of L(3 — 1100) estimated from our 
models are 5-15% lower than the ones for the GLOB- 
ALBKG. This is mainly caused by a decrease in the 
160 /tm flux due to the substantial diffuse emission at 
that wavelength. The diffuse component accounts on av- 
erage for 36% of the total 8 /tm emission of the 160RES 
regions, 23% at 24 /im, 26% at 70 /im, and 34% at 
160 /im. The average difference between our luminosi- 
ties (for the LOCALBKG, i.e., after subtracting the dif- 
fuse emission) and those estimated with the Equation in 
IDale k Helol (|2002f) is 10 ± 3%. In this case, there is 
a clear linear correlation where the difference is higher 
for fainter regions: up to 15% for the faintest source 
at i(3 - 1100) ~ 1O 6 8 L and 5% for the brightest at 
L(3-1100) - 10 s L Q . 

These comparisons suggest that the equations derived 
for entire galaxies with high IR luminosities should be 
corrected by small amounts when applying them for 
sub-galactic regions with moderate SFRs (and for en- 
tire galaxies with low IR luminosities). The reason is 
that the emission of the cold dust is relatively brighter 
for pieces of galaxies or galaxies with moderate SFRs and 
dust contents, and the diffuse emission plays an impor- 
tant role. This diffuse emission is extremely relevant for 
the wavelengths around 160 /im, which contribute the 
most to the TIR luminosity, as we show in the following 
discussion. 

4.2.2. Correlating the three MIPS fluxes to L(8 - 1000) 

Relying on our models of dust emission for the dif- 
ferent 160RES regions in M81, we obtained simple rela- 
tionships (using a singular value decomposition method) 
to obtain L(8 - 1000) and i(3 - 1100) from the fluxes 
at the three MIPS wavelengths (after subtracting the 
stellar component) for sub-galactic regions in the range 
10 7 < L(8 - 1000) < 10 8 L . In all the following Equa- 
tions, L(filter) refers to v c s L v in units of erg s _1 or solar 
units, where v^s is the effective frequency of the filter. 
The relations for the GLOBALBKG case, in solar units, 

1 3 

are : 

L(8 - 1000) = (+1.439 ± 0.415) x L(24) + 
+ (+0.814 + 0.099) x £(70) + (+1.1229 + 0.039) x L(160) 

(2) 

13 All the correlations in this Section have been obtained by 
fitting the data for sub-galactic regions (i.e., not using the data 
point for the entire galaxy) and accounting for the uncertainties. 
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£(3 - 1100) = (+1.832 ± 0.420) x £(24) + 

+ (+0.779 ± 0.094) x £(70) + (+1.176 ± 0.031) x £(160) 

(3) 

For both fits, the data show a scatter of 1% around the 
given relationships. 

The same equations for luminosities estimated after 
subtracting the diffuse emission (LOCALBKG case) are: 

£(8 - 1000) = (+0.202 ± 0.510) x £(24) + 

+ (+0.802 ± 0.160) x £(70) + (+1.303 ± 0.061) x £(160) 

(4) 

£(3 - 1100) = (+0.294 ± 0.551) x £(24) + 

+ (+0.766 ± 0.190) x £(70) + (+1.392 ± 0.063) x £(160) 

(5) 

The scatter of the data points around these relations 
is approximately 4% in both cases. 

4.2.3. Correlating monochromatic MIR fluxes to 
£(8 - 1000) 

Given the poorer resolution of the MIPS 70 and 160 [im 
channels relative to the 24 //m channel, it would be very 
convenient to use the 24 /im emission alone to obtain re- 
liable IR luminosities. This would allow us to study the 
star-forming complexes in M81 (or any other resolved 
galaxy) at high enough angular resolution to resolve in- 
dividual HII regions and compare the IR results with 
other SFR estimators, such as Ha or the UV. 

The degree of reliability of using only one monochro- 
matic flux to estimate the 8-1000 /im luminosity is shown 
in Figure EH including the diffuse emission at all wave- 
lengths (left panel) and removing it (right panel). Here 
we plot the TIR luminosity of each individual region es- 
timated with the three MIPS fluxes (after removing the 
stellar component) as a function of the monochromatic 
luminosity in those MIPS bands. The best fits to the 
GLOBALBKG data are the following (with all luminosi- 
ties in solar units): 

log [L(8 - 1000)] = (+0.63 ± 0.11) + 
+ (+1.020 + 0.015) x log [£(8.0)] (6) 

log [£(8 - 1000)] = (+1.25 ± 0.06)+ 
+ (+0.997 + 0.010) x log [£(24)] (7) 

log [£(8 - 1000)] = (+1.71 ± 0.46)+ 
+ (+0.839 ± 0.065) x log [£(70)] (8) 

log [£(8 - 1000)] = (-0.27 ± 0.36) + 

+ (+1.064 + 0.048) x log [£(160)] (9) 

The diffuse emission has a very significant effect on the 
previous correlations, especially for the 8 and 160 /im 
data. The relationships between the monochromatic 
fluxes and the TIR luminosity for the photometry where 
the diffuse emission has been removed (LOCALBKG 
case) are: 

log [£(8 - 1000)] = (+0.95 ± 0.13) + 



+ (+0.980 + 0.019) x log [£(8.0)] (10) 

log [£(8 - 1000)] = (+1.88 ± 0.27)+ 

+ (+0.897 + 0.041) x log [£(24)] (11) 

log [£(8 - 1000)] = (+1.78 ± 0.27)+ 
+ (+0.824 + 0.038) x log [£(70)] (12) 

log [£(8 - 1000)] = (-0.14 ± 0.30)+ 
+ (+1.050 ± 0.041) x log [£(160)] (13) 

Figure |3| clearly shows that the best estimator of the 
TIR luminosity (the one presenting the lowest scatter) is 
the 160 /tm emission, which accounts for more than 50% 
of the TIR luminosity (see also Equations (2J and 0J|. The 
scatter around a linear relationship for the monochro- 
matic 160 fim data is 5% of the TIR luminosity (6% 
if we remove the diffuse emission). The 24 /tm data 
alone is a much poorer estimator of the TIR luminos- 
ity, presenting a scatter around the given relationship 
of 12% (13%). The goodness of the 70 fim and the 
8 /im data as estimators of the TIR luminosity lies in 
between: 9% (11%) and 8% (7%) scatter around the 
given relationships, respectively. It is important to no- 
tice that, although 8.0 /urn seems to be a better TIR 
luminosity estimator than 24 /im based on the scatters 
of the fits, the IRAC band is more affected by the stellar 
emission contamination (more than a factor of 2 differ- 
ence). Moreover, the 8.0 /im image shows a brighter 
diffuse emission component than the 24 /jm image, 36% 
of the total emission for the former and 23% for t he lat - 
ter (cf. Section |4~2~T|) . According to lFazio et alJ l|2004j) . 
the scattered light in the 8 /im array can only account 
for ^2% of that fraction. This clearly shows that the 
8 (jm emission depends not only on the local star for- 
mation rate but also in some other mechanisms related 
to t he formation and de struction of aromatic molecules 
(see lCalzetti et aLll2005t and references therein) . Indeed, 
in the LOCALBKG case (where the starlight contamina- 
tion and the diffuse emission of uncertain origin should 
be removed), the correlation between the 8 /jm and the 
TIR luminosities improves. 

Figure [21 also shows one interesting point about the 
contribution of the various wavelengths to the integrated 
luminosity. As we move to higher TIR luminosities, the 
contribution of the 160 /im luminosity to the integral be- 
comes less important (as the positive slope shows) , while 
the contribution of the hotter dust emitting at shorter 
wavelengths rises (as shown by the negative slopes of the 
linear fits for 24 and 70 /im). Figure [21 shows that the 
dust contributing the most to £(8 — 1000) has a tem- 
perature of around T ~ 19 K (its emission peaks at 
A - 160 /on) for £(8 - 1000) < 10 8 £ Q . The fraction of 
£(8 - 1000) emitted by the dust at T > 40 K (with emis- 
sion peaking at A < 70 /im) increases as £(8 — 1000) gets 
brighter (i.e., as the current SFR increases). The hotter 
dust contributes as much as the colder dust component at 
£(8 - 1000) > 1O 8 5 £ (i.e., for SFR >0.05 A^yr^ 1 ). 
The fraction of £(8-1000) emitted by dust at T > 125 K 
(with emission peaking at 24 /tm) is always lower than 
10% for £(8 - 1000) < 1O 9 O £ . Another interesting 
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Fig. 3. — Relationship between the MIR and FIR monochromatic emissions and the 8-1000 /an luminosities for star-forming complexes 
in M81 (selected in 160RES resolution images). On the left panel, the data refer to the photometry calculated with a global background 
(GLOBALBKG case in the text). On the right, a local background (LOCALBKG case) was used to remove the diffuse component of the 
emission in the IR. For both panels, the integrated luminosities have been calculated by fitting models of dust emission to the SEDs for 
wavelengths redder than 5 fim. The best linear fits to the data at each wavelength are also shown. The different colors refer to the distinct 
types of regions denned in the left panel of Figure In both panels, average errors for the different luminosities (horizontal bars) and 
luminosity ratios (vertical bars) are plotted to one side of the corresponding data points. 



point in Figure El is that the regions dominated by dif- 
fuse emission (red symbols) seem to deviate significantly 
from the star- forming complexes (blue symbols), show- 
ing a larger contribution to L(8 — 1000) from the colder 
dust. 

Finally, it is also worth noticing that the entire galaxy 
and the nuclear and circumnuclear zones deviate signif- 
icantly from the relationship found for the star-forming 
regions in the arms of M81. These regions present sys- 
tematically hotter dust temperatures than the ones ob- 
tained for the zones in the arms (see Section FOjl . which 
should be linked to other sources of heating, such as the 
central AGN, the bu lge old stars, or dynamical friction 
(IGordon et al.ll200l . 

4.2.4. Correlating MIR colors to L(8 - 1000) 

In Figure QJ we explore the possibility of using a color 
to improve the relationship between the 24 //m and the 
TIR luminosities. In principle, we could expect a color 
involving MIR fluxes to be correlated with the behavior 
of the SED in the FIR, i.e., the amount and emissivity of 
the cold dust (dominating the FIR luminosity) might be 
linked with the properties of the warm dust if both dust 
components share the same heating source. If that were 
the case, we would be able to reduce the scatter in the es- 
timation of L(8 — 1000) from single MIR fluxes. Figure^] 
shows that there is a correlation between the 8.0 ixm-to- 
24 /im color and the TIR luminosity, although the scatter 
is considerable. The best fit to the data (GLOBALBKG 
case) gives: 

log [L{8 - 1000)/L(24)] = (+0.94 ± 0.05) + 
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-0.644 ± 0.102) x log [L(24)/L(8.0)] 



(14) 
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For the data where the diffuse component has been re- 
moved (LOCALBKG case), the correlation is very simi- 
lar: 

log [L{8 - 1000)/£(24)] = (+0.93 ± 0.06)+ 

+ (-0.695 ± 0.155) x log [L(24)/L(8.0)] (15) 

Using this relationship, we are able to reduce the scat- 
ter in the determination of L(8 — 1000) to less than 6% 
(11% for the photometry with no diffuse emission), i.e., 
the estimates are two times better than the ones obtained 
using 24 fim alone. However, the correlation has a sig- 
nificant number of outliers, most of them correspond- 
ing to regions where the star formation is not very high 
and the diffuse emission is very important (inter-arm re- 
gions and zones between the arms and the nucleus). Our 
correl ation lies very near t he predictions from the mod- 
els of [Dale & Hclouj <|2002(1 . and it is significantly differ- 
ent from the relationship found bv lCalzetti et alJ ( 2005) 
for 21 regions in M51, although the slope is very simi- 
lar. These differences must arise from the fact that the 
8-to-24 /mi color depends on the strength of the aro- 
matic emission (contributing mainly at 8 fiva), which 
presents significant variations from galaxy to galaxy de- 
pending on parameters such as the abundance of aro- 
matic molecules, composition (linked to the metallic- 
ity, see lEngelbracht et alJl20 051. or temperature of the 
galaxy fsee lDale et al.M2005|) . Therefore, one should be 
cautious in using Equations 1141 and 1151 for any galaxy 
before it is checked for a larger sample. 

All the correlations found in this Section are summa- 
rized in Table |31 

4.3. Spatially resolved properties of the dust in M81 

Figure [5] shows the MIPS colors as a function of ra- 
dius (de-projected radial distances) for the regions se- 
lected from the 160RES dataset 14 . In our models, the 
radial dependence of the MIPS colors manifests as a 
radial gradient in the dust temperature. The model 
dust temperatures as a function of radius are plotted 
in the insets of Figure This figure shows that the 
IR emission in M81 is dominated by a cold dust com- 
ponent having temperatures around 18 K (the average 
and standard deviation are < T c >= 18 + 2 K). The 
warmer dust component has an average temperature of 
< T w >= 53 ± 7 K. In com parison with these tempera- 
tures, |Devereu3^0i2 l)1995j) obtain an average dust tem- 
perature of T — 30 K (approximately the average of our 
two components) for the entire galaxy and just using the 
IRAS 60 and 10 iim fluxes and mod els with a single 
dust component. iDevereux et alJ l)1995|) also provide ex- 
pected values for the temperature profiles of the dust. 
The temperatures in these profiles range from T ~ 15 K 
to T ~ 50 K (depending on the composition of the dust), 
values that are very similar to the ones obtained for our 
fits. Our dust temperatures are also comparable to those 
estimated from other work based on FIR and sub-mm 
data, which find cold dust temperatur es of 15 — 25 K and 
warm dust temperatures of 40 — 60 K lPirnng_el_aJJ2QQII 
Dunne &: Ealesl l200ll IPopescu et alJ 120021 IBendo et all 
2003HRegan et al.ll2004HEngelbracht et al.ll2004j) 

14 The results in this Subsection refer to the GLOBALBKG case, 
since we are studying the properties of all the dust emission in M81, 
including the diffuse component. 



The cold dust component shows a clear trend of de- 
creasing temperature with increasing radius, although 
the scatter and fitting errors are large. The tempera- 
ture profile is consistent with the observed radial gra- 
dient of the L(160)/L(70) color, which increases from 
L(160)/L(70) - 2 to L(160)/L(70) = 3-4 within the 
disk. For the warm dust component, the temperature 
radial gradient is very small, resulting from a change in 
the L(160)/L(24) color of a factor of 2 from the center to 
the outskirts of M81. The central regions of M81 show 
slightly hotter dust temperatures in the cold component: 
T c = 22.8 ±1.0 K and T c =21.6 + 0.5 for the nuclear and 
circumnuclear regions, respectively. As we mentioned 
earlier, this could be due to the presence of another effi- 
cient source of dust heating, apart from the young stars: 
an AGN for the nucleus and/or older stars (whose den- 
sity is higher in the bulge of the galaxy and may present 
radial changes in the ste llar SEDs) for both zones (see 
IGenzel fe Cesarskvll2000l and references therein). These 
effects could also be responsible for the radial gradient 
in the cold dust temperatures. 

In combination with the radial gradient of the dust 
temperature, our models also reveal a decrease in the 
dust mass surface density as we move to higher radial 
distances (beyond —300"). This is shown in Figure 
where we use de-projected radial distances and surface 
densities. This figure also shows the profile through the 
major axis of M81 for the mass surfac e density of the 
atomic and molecular ga s (estimated by |Djaej£ux ei_al 
1995 using the data in iGottesman fc A^h^hewiQ975l 
and iRots fc Shanel 119751 for the at omic gas, and using 
the 12 CO intensities measu red by iBrouillet et al.1 119911 
and lSage fc Westpfahllll99lT) . For the entire galaxy, the 
dust-to-gas ratio is 0.005 (accounting for both atomic 
and molecular gas). We detect a factor of 10 increase 
in the dust mass surface density from the center of the 
galaxy to the inner part of the arms (at a radial distance 
of approximately 200"). The atomic gas surface density 
also increases by a similar amount in the same region. 
Within the arms (from -200" to -400"), the dust mass 
density is about 10 5 9 A4© kpc -2 (with a scatter of about 
0.2dex). The regions with the highest dust mass surface 
densities are found at radial distances of 300 — 350", just 
where the profiles for the molecular and atomic gas sur- 
face densities peak. The dust mass surface density de- 
creases (by roughly a factor of 10) from —300" to the 
outer regions. This gradient in the dust density should 
produce a radial decrease in the extinction of the stellar 
light, which is confirmed in Section (see Figure ITTfl ■ 

One caveat to our results on the radial gradient of the 
dust mass density arises from the lack of data at wave- 
lengths longer than 160 /urn 15 , where the peak of the 
SED must be located. Our data are not sufficient to ac- 
curately delimit the position of that maximum, which 
results in relatively large uncertainties on the amount of 
dust and emission arising from dust colder than approxi- 
mately 20 K. This means that our estimations of the dust 
mass should be considered as lower limits. More data in 
the sub- millimeter range (at least, at wavelengths red- 
der than —200 fj,m, the range that will be covered by 

15 The only data available for M81 in the sub-mm is a 
JCMT/SCUBA observation of the central 2' of the galaxy which 
do not cover any of the star-forming complexes in the arms. 
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Fig. 5. — Radial plots (using de-projected radial distances) of the MIPS colors for the regions selected at 160RES. The different colors 
refer to the distinct types of regions defined in the left panel of Figure The average surface density of the entire galaxy has been plotted 
at r = 20", and the surface density for the nuclear and the circumnuclear regions has been plotted at r = 40" (slightly shifted from r = 
for clarity). The insets in each panel show the dust temperatures for the cold (left panel) and warm (right panel) components in our models. 
Average errors for the colors and temperatures are shown in each panel. 



TABLE 3 

Summary of correlations between mid-infrared monochromatic luminosities and colors and the integrated 

IR luminosities. 



Estimator 


BKG 


A 


B(8.0) 


B(24) 


B(70) 


B(160) 


a 


L(8 - 1000) 


GLOBAL" 






+1.439±0.415 


+0.814±0.099 


+1.1229±0.039 


1% 


log [L(8 - 1000)] 


GLOBAL 


+0.63±0.11 


+1.020±0.015 








8% 




GLOBAL 


+1.25±0.06 




+0.997±0.010 






12% 




GLOBAL 


+1.71±0.46 






+0.839±0.065 




9% 




GLOBAL 


-0.27±0.36 








+1.064±0.048 


5% 




GLOBAL 


+0.94±0.05 


+0.644±0.102 


+0.356±0.102 






6% 


L(8 - 1000) 


LOCAL" 






+0.202±0.510 


+0.802±0.160 


+1.303±0.061 


4% 


log [L(8 - 1000)] 


LOCAL 


+0.95±0.13 


+0.980±0.019 








7% 


LOCAL 


+1.88±0.27 




+0.897±0.041 






13% 




LOCAL 


+1.78±0.27 






+0.824±0.038 




11% 




LOCAL 


-0.14±0.30 








+1.050±0.041 


6% 




LOCAL 


+0.93±0.06 


+0.695±0.155 


+0.305±0.155 






11% 


L(3 - 1100) 


GLOBAL" 






+1.832±0.420 


+0.779±0.094 


+1.176±0.031 


1% 


L(3 - 1100) 


LOCAL" 






+0.294±0.551 


+0.766±0.190 


+1.392±0.063 


4% 



Note. — The Table shows the coefficients of the correlations (in logarithmic scale except for the correlations in- 
volving the 3 MIPS wavelengths) between the MIR and FIR monochromatic luminosities and the total IR luminosity 
integrated between 8 and 1000 fim [1/(8 — 1000)] and 3 and 1100 (im [L(3 — 1100)]. The equation for each correlation 
is: log [L(8 — 1000)] = A + B(Aj) X logL(Ai) [except for the correlation between the 3 MIPS luminosities and the 
integrated IR emission, where the correlation is L(8 — 1000) = B(Ai) X L(Ai)]. The sum extends over the wavelenghts 
of the last IRAC channel and the 3 MIPS bands (8.0, 24, 70, and 160 /im). In the last column, a represents the scatter 
of the data around the given relation. All luminosities are given in erg s _1 . 
"This correlation is given in linear scale, not logarithmic. 



Herschel) are necessary to account for the coldest dust 
component and robustly estimate the total amount of 
dust. 

5. COMPARISON OF STAR FORMATION RATE 
ESTIMATORS 

In this Section, we compare different SFR tracers for 
59 individual HII regions in M81 with available extinc- 
tion measurements. All of these HII regions are resolved 
in our 24RES image set (see the right panel of Figure 



The sample of 59 sources includes the 18 HII regions 
in the disk and outskirts of M81 f at radial distances 
betwe en 3 and 15 kpc) for which iGarnett fc Shields! 
(1987) obtained optical spectra. These authors mea- 
sured Balmer decrements that, appropriately corrected 
for stellar absorption and temperature effects, provide a 
good estimate of the extinction of the gas emission. The 
remaining 42 sources in our sample are giant HII regions 
(at radial distances between 3 and 10 kpc) for which 
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Fig. 6. — Radial profile (using de-projected radial distances) of 
the dust mass surface density (also corrected for inclination). The 
different colors refer to the distinct types of regions denned in the 
left panel of Figure The average surface density of the entire 
galaxy has been plotted at r = 20", and the surface density for the 
nuclear and the circumnuclear regions has been plotted at r = 40" 
(slightly shifted from r = for clarity). The profiles through the 
major axis for the mass surface density of t he atomic gas (contin- 
uous line) and molecular gas (dashed line, Dcvcrcux et al.1119951) 
are also plotted, with the scale shown on the right vertical axis. 



iKaufman et alJ l|1987fl obtained Ha and radio (at 6 and 
20 cm) observations, estimating optical extinctions for 
most of them. These extinction estimations may be af- 
fected by radio synchrotron emission not d irectly linked 
to the recent star formation (fCondo n ll 19921) . 

In the following, we compare the luminosities of the 
HII regions in M81 for 5 SFR tracers: the UV emis- 
sion, the Ha nebular emission, and 3 luminosities in the 
IR wavelength range [L(8 - 1000), £(24), and L(8), all 
referring to the photometry measured with a LOCAL- 
BKG]. In this comparison of SFR estimator s , we also 
include the results obtained bv lCalzetti et al.1 <)2005D for 
the central region of M51. In that work, Ha and Paa 
narrow-band fluxes, jointly with IRAC and MIPS emis- 
sions, are measured for HII regions with high extinctions 
(1 < A(V) < 4). The combination of the data for M81 
and M51 allows us to probe a wide range of attenuations 
(from A(V) - to A(V) ~ 4) and SFRs (from 3 x 10~ 4 
to 0.1 A^oyr" 1 , or 5 x 10~ 3 to 0.3 M Q yr^kpc -2 ). 

The extinction corrected Ha luminosities used in this 
Section have been estimated by adopting the following 
differential values between the Ha, H/3, and Paa at- 
tenuation factors, fc(A) 16 : fc(H/3) — fc(Ha) = 1.23 and 
fc(Ha) - fc(Paa) = 2.02, with fc(Ha) = 2.41. These 
factors correspond to a screen attenuation law with 
A(V) ~ 0.5 for M81 and A{V) ~ 3.0 for M51, homo- 

16 The attenuation at a given wavelength, A(X), is given by 
A(A) = fc(A) x E(B - V) 



geneous distribution of the dust, and dust properties 
typic al of the Small Magellanic Cloud ijWitt fe; Gordonl 
2000). These attenuation laws see m to best repro duce 
the obse rvations for M81 (se e, e.g.. lHill et al.lll995|) and 
M51 (see JCalzettiet^L| l2005|) . Other attenuation laws in 
IWitt & Gordonl <|2000|) with different star-dust geome- 
tries a nd distributions, or extinction laws such as the 
ones in lCardelli et alJ (^989) , provide similar factors (dif- 
ferent by less than 1%). We have assumed the theoret- 
ical line ratios Ha/H/3 =2.86 and Ha/Paa=8.46, values 
expected for a low density gas (n c ~ 10 2 cm -3 ) with 
T c = 10 4 K in the recombinatio n Case B including colli- 
sional effects (Osterbrock 1989). 

5.1. The TIR luminosity as a SFR tracer 

The left panel of Figure [7| shows the comparison be- 
tween the TIR luminosity (calculated with Eauation ll5fl 
and the extinction corrected Ha emission for the HII re- 
gions in M81 and M51. The best linear fit to the data 
for both galaxies is 17 : 

log{£(8 - 1000)} = 

(+2.3 ± 1.2) + (+0.992 ± 0.031) x log{i corr (ffa)} (16) 

The scatter of the data around this relationship is 30%. 
Th is equation is very d ifferent from the one obtained 
bv iCalzetti et all (2005) for the M51 regions alone. In 
fact, Figure [7| (see the inset in the left panel) shows that 
the M81 and M51 data points seem to follow different 
trends. The slopes obtained for the M81 or the M51 
data alone are very similar (~ 0.9), but there is an offset 
in the zero point (shown in the inset by the offset in the 
cloud of points for both galaxies) . For the same value of 
L corr (i?a), the M51 regions are brighter than the M81 
ones in the IR. This offset is consistent with the larger 
A(V) values found for the M51 central region in com- 
parison with the extinctions found for the HII regions in 
the arms of M81. For the regions with larger extinctions 
(the ones in M51), more energy is absorbed by the dust 
and reemitted in the IR. 

Classically, total SFRs for galaxies (or pieces of galax- 
ies) have been obtained in two different ways. One op- 
tion is using an estimator affected by dust attenuation, 
such as an emission-line or the UV emission, and cor- 
recting the observed flux for the effects of extinction. 
Another possibility is directly using a SFR tracer that is 
not affected by dust attenuation, such as the IR or radio 
emission. Both ways of estimating total SFRs were com- 
pared in the left panel of Figure \7\ (where we compare 
extinction-corrected Ha luminosities with TIR luminosi- 
ties). In principle, the IR emission is linked to the dust 
heated by the light coming from the newly-formed stars 
in the HII regions. The UV and Ha emissions (without 
any extinction correction) are linked to the photons aris- 
ing from those young stars that do not interact with the 
dust at all. This means that we can obtain an estimate 
of the SFR from the observed Ha or the UV emission 
that cannot be detected by the IR because these pho- 
tons do not interact with the dust, and on the other 
hand, we can account for the emission-line or UV pho- 
tons which are extinguished by dust by using IR obser- 

17 All the linear correlations in this Section will be given for 
luminosities in units of ergs -1 . 



15 




log(L corr (Ha) [erg s -1 ]} log{L Qbs (H«) [erg s -1 ]} 

Fig. 7. — On the left, comparison between the extinction corrected Ha and the TIR luminosities for selected HII re gions in M81 (with 
differ ent colors a nd symbols referring t o the distinct types of regions defined in the right panel of Figure^ extracted from lGarnett fc' Shiclds 
119871 -GS87- and lKaufman et aH1987l -K87-) and in M51 (black triangles; data extracted from Calz ettiet alJ2005ll . On the right, the TIR 
luminosity of the same sample of HII regions is compared with the Ha luminosity absorbed by the dust. For both panels, we show the 
typical uncertainties linked to each catalog of HII regions in M81. The best linear fit is plotted with a continuous gray line. The insets 
show the same plots but using a luminosity ratio in the vertical axis to enhance the differences in the trends followed by the datasets for 
each galaxy. The dotted line in the right panel (and the ins et) shows the lum inosity relation for slope unity and a ratio between the Ha 
and TIR luminosities given by the SFR calibrations found in lKennicutU 0998). 



vations. Consequently, an alternative way to obtain a to- 
tal SFR would be to add the 'unobscured SFR' derived 
in the UV/optical (without any extinction correction) 
and the 'obscured SFR' obtai ned from the IR luminos- 
ity (see iKennicutt et al J 1200^1 . This total SFR should 
be the same as the SFR derived with the UV/optical 
estimators after correcting them for the effects of extinc- 
tion and taking into account the spatial variations of the 
dust content in star-forming regions, which might result 
in the UV/optical and IR emissions arising from different 
regions (with different attenuations). It should also coin- 
cide with the SFR derived with the IR emission alone, if 
we can obtain a calibration for it that takes into account 
the amount of star formation not obscured by dust (or if 
this star formation is negligible in comparison with the 
star formation hidden by dust). 

The alternative method of obtaining the total SFR 
adding the estimations of the unobscured SFR (obtained 
with the observed Ha or/and UV luminosities) and the 
obscured SFR (obtained with the IR emission) requires 
the existence of a (simple) correlation (close to pure pro- 
portionality) between the extincted correction for the 
UV/optical SFR tracers and the emission in the IR. The 
right panel of Figure [7| explores this correlation for the 
TIR luminosity, comparing it with the Ha luminosity ex- 
tinguished by dust (L abs (Ha), calculated by subtracting 
the observed luminosity from the extinction corrected lu- 
minosity) for the HII regions in M81 and M51. The best 
linear fit to the data is: 

log{£(8 - 1000)} = 

(+8.9 ± 1.0) + (+0.826 ± 0.025) x log{L ahs (Ha)} (17) 

The scatter of the data around this relationship is 30%. 
The scatter of the points for the M81 data is larger than 



for the M51 data (45% and 20%, respectively). Most of 
this scatter (and the scatter observed in all the panels in 
Figures [7| through |5J is due to the larger uncertainties in 
the extinctions derived for the HII regions in M81 in com- 
parison with the Paa-to-Ha derived extinctions for M51. 
The scatter might also be caused by spatial changes in 
the extinction properties within an individual HII region 
(i.e., the actual net extinction and the attenuation law). 
If the central parts of the HII regions had larger extinc- 
tions than the outskirts, our values of L corr (Ha) would 
be overestimated. To test this effect, we measured fluxes 
for the same HII regions but with smaller apertures. The 
correlation was slightly better in this case, with a scatter 
around the linear relationship of 25%. 

The different trends shown by the M81 and M51 data 
in the inset of the left panel of Figure [7| disappear in 
this plot, although the scatter of points is slightly larger, 
mostly because of the M81 points (and especially the 
ones with extinctions estimated with radio data) . In ad- 
dition, the slope of the correlation differs significantly 
from unity (0.826+0.025), the value that is expected if 
the two quantities are closely correlated. This behavior 
suggests that the heating of the colder dust, which dom- 
inates the TIR luminosity, is not directly caused (ex- 
clusively) by the Lyman photons. The right panel of 
Figure {7\ also shows the luminosity relation with slope 
unity and a ratio between the Ha and TI R luminosi- 
ties gi ven by the SFR calibrations found in iKennicuttl 
(1998). Note that Kennicutt's calibration was built for 
entire starburst galaxies assuming continuous star for- 
mation for the last 10-100 Myr. Despite this, the rela- 
tion fits the data points for the M81 HII regions very 
well, but it departs significantly from the brightest M51 
points. This effect indicates differences in the attenua- 
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tion properties (for example, in the attenuation law) of 
the mildly extincted M81 regions in comparison with the 
heavily attenuated M51 central zone. 

5.2. The rest-frame 24 fJ-m luminosity as a SFR tracer 

For the nuclear region of M51, iCalzetti et all l)2005fl 
found that the 24 /im emission correlates tightly (even 
better than the TIR luminosity) with the extinction cor- 
rected Paa luminosity, suggesting that the warm dust 
responsible for the 24 /im emission is heated by the ioniz- 
ing photons. Indeed, it seems reasonable that the hottest 
dust in a galaxy could be more efficiently heated by 
the very energetic photons arising from the hottest (and 
youngest) stars than by redder photons from older stars. 
In contrast, the colder dust could be heated not only by 
the ionizing flux, but also by the emission arising from 
cooler (and older) stars. In such case, the luminosities 
in the MIR would be a better tracer of the most recent 
star formation than the TIR luminosity, which is dom- 
inated by the emission of colder dust (see Equations 
and This is supported by the fact that the morpholo- 
gies of M81 and M51 at 24 /zm rese mbl e those of the 
Ha emission (see iGordon et al.ll2004l and ICalzetti et aT 
20051 see also lHinz et al.ll2004l for M33, and lHelou eTaT 
20041 for NGC300), which closely traces the location of 



the ionizing stars (in the HII regions). In contrast, the 
morphologies in the MIR are not as similar to the UV 
images, given that the UV flux is largely affected by ex- 
tinction and the emission is normally linked to the star 
formation on a longer timescale. 

We explore the possibility of using the 24 /im emis- 
sion alone as a SFR tracer in Figure |S1 The left panel 
shows the correlation between £(24) and L corr (Ha) for 
the combined sample of HII regions in M81 and M51. 
The best fit to the data is: 

log{L(24)} = 

(-10.7 ±1.4) + (+1.302 ±0.036) x log{L corr (Ha)} (18) 

The global scatter around this relationship is 40%. As 
seen in the left panel of Figure the data points for 
M81 and M51 seem to fall on two different linear rela- 
tions with very similar slopes (very close to unity) but 
offset (see in the inset in the left panel of Figure 
This effect is c onsistent with the d ifferent L(24)/ SFR 
ratios found bv ICalzetti et all (2005) for the HII regions 
in M51 and for entire galaxies, including star bursts and 
ultraluminous infrared galaxies. While the galaxies with 
very intense star formation show larger L(2A)/SFR ra- 
tios than the central region of M51, the HII regions in 
M81 show smaller ratios (on average, a L(24) / SFR ratio 
a factor of three lower than the centr al region of M51). 
As pointed out by Calz etti et alJ l|2005() . £(24) alone does 
not seem to be a universal estimat or of the total SFR of 
a gala xy or HII region. However, lAlonso-Herrero et al.l 
(2006) have derived an empirical relation that provides 
an accurate fit to the SFR (obtained from Paa data) to 
24 /im luminosity for the dusty M51 HII regions upward 
in the luminosity through L(8 - 1000) ~ 10 12 L@. 

In the right panel of Figure |H1 we show the alternative 
calibration of the 24 /im luminosity as a tracer of the 
amount of Ha photons extinguished by dust, i.e., not of 
the total star formation, but only of the star formation 
which is undetectable by the observed Ha emission due 
to dust extinction. The best fit to the data is: 



log{Z(24)} = 

(+1.2 ± 1.3) + (+1.002 ± 0.034) x log{L abs (Ha)} (19) 

The global scatter around this relationship is 35%. The 
most deviant points correspond again to the HII regions 
in M81 with a radio-based extinction. The slope of the 
fit is consistent with a value of unity, which favors a 
direct link between the Lyman photons (or the emission 
arising from the stars that dominate the production of 
Lyman photons) and the heating of the hottest dust in 
M81 and M51. Note that the SFR calculated in this way 
should be added to the SFR estimated from the observed 
emission-line to obtain the total SFR. 

5.3. The 8 /im luminosity as a SFR tracer 

Figure [5] shows the comparison of the 8 /im luminosity 
(corrected for the stellar light contamination) with the 
extinction corrected Ha luminosity and the Ha luminos- 
ity extinguished by dust. To be able to compare with 
our data, the 8 /im emission for the regions in M51 has 
been corrected for the effect of the calibration of extended 
sources in IRAC images (a correction of 0.133 dex). The 
best linear fits to the data are: 



(-4.5 ± 1.6) 



log{L(8)} = 
-1.147 ±0.042) x log^^XHa)} (20) 

log{L(8)} = 



(+5.2 ± 1.2) + (+0.902 ± 0.031) x log{L abs (Ha)} (21) 

The scatter around these relations is larger than 60% 
in both cases. Comparing with Figures and |H1 Figure^l 
shows that the 8 /im emission is a poorer estimator of the 
total or extincted SFR than the 24 /im or TIR emissions. 
Within each galaxy, the data seem to follow a linear re- 
lation (with a slope that differs significantly from unity) , 
but there are strong differences between the two galaxies. 
Three factors can be responsible for these differences: 1) 
an important residual contamination in the 8 /im fluxes 
from stellar emission, which might depend on the size of 
the apertures used and the properties of the stellar pop- 
ulation enclosed; 2) galaxy-to-galaxy variations in the 
metallicity and aromatic properties and spectrum; and 
3) a heterogeneous nature in the agents responsible for 
the excitation of the arom atic molecules that d ominate 
the 8 /im emission (see also ICalzetti et al"ll2005l and ref- 
erences therein). 

All the correlations found in this Section are summa- 
rized in Table El 

5.4. Summary of the IR SFR tracers 

All the results for M81 and M51 described in the pre- 
vious subsections support the idea that the heating of 
the dust emitting at 24 /im (the dust at T = 40 - 50 K) 
is dominated by the light arising from the same stars 
that produce most of the ionizing flux, which is directly 
correlated with the emission-line fluxes. In addition, the 
spectrum around 24 /im shows a remarkable absence of 
any bright emission or absorption feature (such as aro- 
matic molec ules) f or a wide var iety o f galaxies (see, e.g., 
iSturm et al.ll2000l lArmus et al1l200l iSpoon et al.ll200l 
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logfL corr (Ha) [erg s -1 ]} log{ L abs (Ha) [erg s -1 ]} 

Fig. 8. — On the left, comparison between the extinction corrected Ha and the 24 (im luminosities for selected HII regions in M81 (with 
differ ent colors a nd symbols referring t o the distinct types of regions defined in the right panel of Figure^ e xtract ed from|Garnctt & Shields! 
119871 -GS87- and lKaufman et alJU987l -K87-) and M51 (black triangles; data extracted from lCalzetti et aljl20031 . On the right, the 24 fim 
luminosity of the same sample of HII regions is compared with the He* luminosity absorbed by the dust. For both panels, we show the 
typical uncertainties linked to each catalog of HII regions in M81. The best linear fit is plotted with a continuous gray line. The insets 
show the same plots but using a luminosity ratio in the vertical axis to enhance the differences in the trends followed by the datasets for 
each galaxy. The dotted line in the right panel (and the inset) shows the luminosity relation for slope unity. 




log{L corr (Ho() [erg s -1 ]} log{L abs (Ha) [erg s" 1 ]) 

Fig. 9. — On the left, comparison between the extinction corrected Ha and the 8 /an luminosities for selected HII re gions in M81 (with 
differ ent colors a nd symbols referring to the distinct types of regions defined in the right panel of FigureHl extracte d from]Garnctt & Shields! 
119871 -GS87- and lKaufman et al.M1987l -K87-) and M51 (black triangles; data extracted from lCalzetti et al.ll2005l) . On the right, the 8 fira 
luminosity of the same sample of HII regions is compared with the Ha luminosity absorbed by the dust. For both panels, we show the 
typical uncertainties linked to each catalog of HII regions in M81. The best linear fit is plotted with a continuous gray line. The insets 
show the same plots but using a luminosity ratio in the vertical axis to enhance the differences in the trends followed by the datasets for 
each galaxy. 



ISmith et all 120041: iHouck et all 12004). This means that 
the 24 /Ltm band is adequate to directly measure the con- 
tinuous emission from warm dust (heated by ionizing 
photons) without line contaminants, in contrast to the 
significant contamination of other parts of the IR spec- 
trum by aromatic emissions (e.g., 8 or 12 /im) or/and 



emission from dust heated by sources not directly linked 
with the recent star formation (e.g., 160 fim). All these 
reasons support the idea that the 24 /im emission alone is 
a good estimator of the obscured SFR for HII regions in 
a wide range of extinctions (from A(V) ~ to A(V) ~ 4) 
and SFR (from 10~ 4 to 0.3 A^gyr 1 ). In the same 
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sense, the 8 /im emission is a poorer estimator of the 
SFR. But the IR luminosities alone are not directly pro- 
portional to the total SFR for HII regions in different 
galaxies. Regions with higher extinctions present la rger 
SFRs CseelHopkins et aDl20Qll Isiillivan et al.ll20Qll and 
iPerez-Gonzalez et all2003dj) and higher IR-to-Ha ratios 
that increase non- linearly with L c °"(Ha). There are two 
reasons for this behavior: 1) the absorption of photons 
by dust does not scale linearly with the radiation field 
seen by the dust; 2) as we increase the extinction, larger 
zones of the HII regions become optically thick and the 
emission-line observations are not able to recover all the 
emitted flux. More data for other galaxies (spanning 
more extinction, SFRs, metallicities, and stellar popula- 
tion properties) should be added to this study to prove 
the universality of the correlations given in the previous 
subsections. 

5.5. Calibration of the UV luminosity as a SFR tracer 
compatible with Ha and the IR 

In this Subsection, we compare the UV-derived SFRs 
with the Ha observations and obtain a calibration of the 
UV data as a SFR tracer compatible with the Ha emis- 
sion. This calibration will allow the use of UV, Ha, and 
IR observations to analyze in detail the SFRs and ex- 
tinction properties of HII regions. 

The relationship between the luminosity in the FUV 
or NUV and the SFR is well known, given that the UV 
emission arises mainly from hot massive stars. However, 
the emission at these wavelengths is strongly affected by 
dust attenuation. This attenuation has been commonly 
estimated using the slope in the UV part of the spectrum. 
The UV slope was calibrated against the actual extinc- 
tion in the UV, measured with the ratio between the 
IR and UV emiss i ons, f or a sample of starburst galaxies 
bv iMeurer et alJ l[l999). This calibration has been ex- 
tensively used in the literature, but several studies have 
shown that the relationship found for starbursts cannot 
be applied to all g alaxies (see, e.g.. lBell2002tlKong et alJ 
12004 iBuat et all 120 05) . Star -forming regions in M81 
(|Gordon et al.l2004D and M51 ijCalzetti et all2005ft have 
also been shown to depart significantly from the star- 
burst relationship. All these works have shown that the 
relationship between the UV slope and the extinction de- 
pends on factors such as the age of the stellar population, 
the metallicity and the attenuation law (which depends 
on the dust properties and the dust-stars geometry) . The 
difficulties in extinction correction prevent the use of the 
UV emission alone as an accurate SFR tracer. 

Another possibility is to use the UV emission to es- 
timate the SFR not obscured by dust, and to add this 
component to the obscured SFR estimated with th e IR 
emissi on. This method is similar to the one used in[BeIJ 
( 2003) . It is also analogous to the procedure we followed 
in the previous subsections, where we showed the good 
correlation between the 24 /zm luminosity and the SFR 
that is obscured by dust and cannot be detected with 
an emission-line such as Ha or Paa. The narrow-band 
data for those lines are very expensive to obtain in terms 
of observing time. Therefore, it would be convenient to 
have an estimate of the SFR not hidden by dust using 
other observations which were relatively easier to obtain, 
such as the GALEX UV data. The goal is to translate 
the broad-band UV emission to an integrated emission in 



the UV, and that to a non-obscured SFR consistent with 
the ionized hydrogen emission. In principle, the Ha emis- 
sion traces a stellar population with a shorter timescale 
than the UV (about 10 Myr for Ha and 100 Myr for the 
UV) , given that the emission-line is linked to the hottest 
(and youngest) stars dominating the production of ioniz- 
ing photons. The Ha/UV ratio should be more and more 
independent of the star formation history of the galaxy 
as we move to sh orter wavelengths in the UV ijBuat et alJ 
l2002HBell2002|) . The extinction also plays an important 
role, since it affects the stellar and ionized gas emission 
differently fsee. e.g.. ICalzetti et alJl2?KXi iGharlot fc Fall 
2000; Bc M~eTal1l2002h . 

In Figure we plot the relationships between the 
observed Ha emission and the luminosities in the two 
GALEX UV bands for the HII regions in M81 and the 
central part of M51. The left panel shows that the NUV 
luminosity is not directly proportional to the Ha emis- 
sion (a proportionality would show in the figure as a 
cloud of points around an average value, but there is 
an offset between the M51 and M81 points). The lack 
of correlation is related to the fact that the NUV emis- 
sion arises from cooler stars than the ones that dominate 
the production of ionizing photons (i.e., the difference in 
the timescales of the star formation traced by Ha and 
the NUV is relatively large). The different luminosity 
ratio found for M81 and M51 should also be related to 
the larger (extinction corrected) SFRs observed for the 
regions in M51. The correlation between the FUV and 
Ha emission is considerably better (i.e., they are roughly 
proportional), although there still seems to be an offset 
between the M81 and M51 data. This panel shows that 
the FUV luminosity probes a stellar population that is 
closer to that producing the Lyman photons, as expected. 

The first two panels of Figure EI] show that there is not 
a direct correlation between the broad-band UV emission 
and the observed Ha luminosity, since parameters such as 
the age of the stellar population and the extinction affect 
both quantities in differing amounts. In the last panel of 
Figure 1101 we explore the possibility of using the two 
GALEX luminosities to obtain an estimation of the un- 
obscured SFR consistent with the observed Ha emission. 
For this purpose, we obtained a relationship between the 
observed Ha luminosity and the emission in the two UV 
bands using a singular value decomposition method. By 
using the two UV fluxes, the relation is intended to cope 
with the two parameters affecting the UV slope, the age 
and the extinction (for a fixed metallicity), and to allow 
obtaining an estimate of the unobscured SFR from UV 
data consistent with the observed Ha flux. The result is 
the following: 

log{L(Ha)} = 

+ lA37xlog{L(FUV)}+(-0A78)x\og{L(NUV)} = log{L(GALEX' 

This empirical relationship allows estimating the un- 
obscured SFR using the GALEX data with a scatter of 
40%. More data for other galaxies should be added to 
this study to prove the universality and usefulness of such 
a correlation. 

5.6. The radial trend in the dust attenuation 

Figure ^2 shows a radial plot of the ratio between the 
obscured SFR (calculated with the 24 /urn luminosity, 
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Fig. 10. — In the first two panels, we show the comparison between the observed luminosities in the two GALEX UV bands and the 
observed Ha emission. In the right panel, we com par e the observed Ha emission with a combined UV luminosity obtained with the two 
GALEX bands, L(GALEX) (defined in Equation 1221 . built to match the observed Ha emission. In this panel, the expected luminosity 
ratio of unity (if L(GALEX) and L(Ha) match perfectly) is also shown. 
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Fig. 11. — Radial dependence (de-projected) of the ratio between 
the obscured SFR esti mate d from the 24 fim luminosity (using the 
calibration in Equation ll9l and the total SFR calculated with the 
observed Ha and 24 fim emissions. The different colors refer to the 
distinct types of regions defined in Figure Q The average surface 
density of the entire galaxy has been plotted at r = 20", and the 
surface density for the nuclear and the circumnuclear regions has 
been plotted at r = 40" (slightly shifted from r = for clarity). 



ratio, i.e., there is a decrease in the extinction of the stel- 
lar light as we move away from the center of the galaxy. 
This gradient is consistent with the results we found in 
Section^ about the radial trend in two other quantities: 
the dust temperature and the dust mass surface den- 
sity. It is also consistent with resul ts obtained for M81 
and other galaxies in the literature l| Boissier et al.ll2004t 
iHolwerda et alJl2005t iRoussel et al1l2005HPrescott et alJ 



Figure ^2 also shows that for the entire M81 galaxy, 
the IR just sees ~50% of the total SFR; the other 50% is 
linked with star formation that is not obscured by dust. 
The largest percentages are observed for the regions dom- 
inated by very diffuse emission, where about 65% of the 
star formation is seen in the IR. For the most remote re- 
gions in the disk, the IR alone can miss up to ~80% of the 
star formation. The gradient in the SFR ratio implies a 
radial change in the Ha extinction from A(fta) ~ 1.0 to 
^4(Ha) ~ 0.2 along the disk of M81, the observed range 
for the 24RES regions used in the previous subsections. 

6. SUMMARY AND CONCLUSIONS 

The recent star formation in the early-type spiral 
galaxy M81 has been analyzed in two dimensions us- 
ing imaging observations spanning the electromagnetic 
spectrum from the far ultraviolet (UV) to the far in- 
frared (IR). The data have been compared to models of 
stellar, gas, and dust emission, obtaining results for dif- 
ferent sub-galactic regions whose sizes range from those 
typical of individual HII regions (0.1 kpc) to the sizes of 
dynamical structures (2 kpc). 

The main results of this work are: 



translated to an equivalent absorbed Ha luminosity with 
Equation [HO and from that to a SFR using the calibra- 
tion in iKennicuttll 1 998h and the total SFR (calculated as 
the sum of the unobscured and obscured SFRs with the 
observed Ha and 24 /im emissions) for the regions de- 
fined in the 160RES image set 18 and using the LOCAL- 
BKG photometry. There is a clear gradient in the SFR 

18 We have used the 160RES regions because they span a wider 
range of radial distances than the 24RES regions 



1. - Using stellar population models, we confirm that 

the luminosities in the mid-IR bands (from 3 to 
8 iim) contain substantial contributions from both 
dust and stellar light. These contributions vary sig- 
nificantly depending on the spatial resolution and 
environment, suggesting the existence of mid- and 
far-IR diffuse dust emission not directly linked to 
the recent star formation. 

2. - The IR emission of M81 can be modeled with three 

components: 1) a cold dust component of average 
temperature T c — 18 ± 2 K; 2) a warm compo- 
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TABLE 4 

Summary of correlations between different SFR estimators. 



Estimator 


log [L corr (Ha)] 




log [L abs (Ha)] 






a b 


<j 


a b 


a 


log [L(8 - 1000)] 
log [1,(24)] 
log [L(8)] 


2.3±1.2 0.992±0.031 
-10.7±1.4 1.302±0.036 
-4.5±1.6 1.147±0.042 


30% 
40% 
60% 


8.9±1.0 0.826±0.025 
1.2±1.3 1.002±0.034 
5.2±1.2 0.902±0.031 


30% 
30% 
60% 



Note. — The Table shows the coefficients of the correlations between the luminosity 
given in the first column and log L COTI (Ho) or logL abs (Ha). The equation for the 
correlation is: log [L(A)] = a + b X log [L 1 (He?)]. In the fourth and seventh columns 
of the Table, a is the scatter around the linear relation. All luminosities are given in 
erg s _1 . 



nent with T w — 53 ± 7 K; and 3) polycyclic aro- 
matic molecules with different properties. While 
the emission of the hotter dust (peaking at around 
24 /ttm) is very concentrated in the HII regions, 
the emission of the colder dust component (peak- 
ing near 160 /im) and the aromatic molecules (con- 
tributing importantly to the 8 /im emission) show 
a more diffuse distribution. This suggests the pres- 
ence of an alternative heating source different than 
the stars recently formed in the HII regions (older 
stars in the disk), or the importance of the trans- 
fer of energy (e.g., the ionizing flux) from the HII 
regions to the interstellar medium (i.e., the escape 
of photons from the HII regions that interact with 
dust in distant zones of the disk). 

3. - A negative radial gradient is observed in the 

dust temperature of the cold dust component, the 
dust mass surface density (similar to the gradi- 
ent of_tiie_jnokecular__gas surface density found 
bv lDevereux et aTlll995|) . and the comparison be- 
tween the star formation rate (SFR) indicators af- 
fected by dust extinction (U V and Ha emission) 
and the IR emission (see also lRoussel et al.1l2005t 
lPrescott"et^Ill200r^ . 

4. - For the entire M81 galaxy, about 50% of the star 

formation is obscured by dust and can be seen in 
the IR, and the other half can be traced directly 
with the observed Ha emission. The percentage of 
obscured star formation ranges from 60% in the in- 
ner regions of the galaxy to 30% in the most distant 
zones. 

5. - Based on our dust emission models, we have built 

several empirical relations between the monochro- 
matic fluxes in the IR and the total IR emission [in- 
tegrated between 8 and 1000 /im, L(8 — 1000)], use- 
ful for sub-galactic regions in spiral galaxies with 
10 7 < L(8 - 1000) < 10 8 L Q . These correlations 
differ significantly from the ones found in the liter- 
ature for entire IR-bright galaxies, mainly because 
of the presence of a colder dust component. 

6. - We have compared five different SFR estimators for 

individual HII regions in M81 and M51: the UV, 
the Ha emission, and three estimators based on 
Spitzer IR data. Several empirical relations have 



been built to estimate the total SFR of HII regions 
from IR data alone, or from a combination of IR 
fluxes (which account for the obscured star forma- 
tion) and UV/optical data (directly measuring the 
unobscured star formation). The typical scatter of 
data around these relations is 30-50%, with possi- 
ble important additional variations of the correla- 
tions from galaxy to galaxy. 

7. - Among the SFR estimators in the IR, the 24 //m 

emission has the best correlation with the Ha lu- 
minosity extinguished by dust. The correlation is 
worse for the total IR luminosity, which suggests 
that the colder dust dominating the total IR lu- 
minosity is (also) heated by redder photons. The 
8 fim luminosity is not directly correlated with 
the ionized gas emission, making it in an unre- 
liable SFR estimator. These results are co nsis- 
tent with those found in iCalzetti et al.1 l)2005(l and 
iKennicutt et alJ (|2006f> . 

8. - Important variations from galaxy to galaxy are 

found when estimating the total SFR with the 
24 /mr or the total IR emission alone. We sug- 
gest the combination of the Ha (or the UV) and 
an IR luminosity (espe cially the rest-frame 24 /im 
emission, calibrated in IKennicutt et alJ 120061 and 
in this paper with the M81 and M51 data) to ob- 
tain the most reliable estimates of the total SFRs 
for HII regions (and, by extrapolation, for galax- 
ies). Those emissions probe the unobscured and 
obscured star formation, respectively. 

9. - Additional data for other galaxies (spanning more 

extinction, SFRs, metallicities, and stellar popula- 
tion properties) should be added to this study to 
prove the universality of the correlations found in 
this paper. 



We thank an anonymous referee for her/his useful com- 
ments. Support for this work was provided by NASA 
through Contract no. 1255094 issued by JPL/Caltech. 
This work is part of SINGS, the Spitzer Infrared Nearby 
Galaxies Survey, one of the Spitzer Space Telescope 
Legacy Science Programs, and was supported by the 
JPL, Caltech, contract 1224667. This work is based in 



21 



part on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Labora- 
tory, Caltcch under NASA contract 1407. GALEX is a 
NASA Small Explorer launched in 2003 April. We grate- 
fully acknowledge NASA's support for construction, op- 
eration, and scientific analysis of the GALEX mission. 
This research has made use of the NASA/IPAC Extra- 



galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. P.G. P.-G. and A. G. de P. also 
wish to acknowledge support from the Spanish Programa 
Nacional de Astronorma y Astroffsica under grant AYA 
2004-01676. 



REFERENCES 



Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, 

ApJ, 290, L25 
Alonso-Herrero et al. 2006, ApJ (in press) 
Armus, L. et al. 2004, ApJS, 154, 178 
Bell, E. F. 2002, ApJ, 577, 150 
— . 2003, ApJ, 586, 794 

Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 

Bell, E. F., Gordon, K. D., Kennicutt, R. C, & Zaritsky, D. 2002, 

ApJ, 565, 994 
Bcndo, G. J. et al. 2003, AJ, 125, 2361 

Boissier, S., Boselli, A., Buat, V., Donas, J., & Milliard, B. 2004, 
A&A, 424, 465 

Boselli, A., Lequeux, J., & Gavazzi, G. 2004, A&A, 428, 409 
Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, 

Kauffmann, G., Heckman, T., & Brinkmann, J. 2004, MNRAS, 

351, 1151 

Brocklehurst, M. 1971, MNRAS, 153, 471 

Brouillet, N., Baudry, A., Combes, F., Kaufman, M., & Bash, F. 

1991, A&A, 242, 35 
Buat, V., Boselli, A., Gavazzi, G., & Bonfanti, C. 2002, A&A, 383, 

801 

Buat, V. ct al. 2005, ApJ, 619, L51 
Calzetti, D. 2001, PASP, 113, 1449 

Calzctti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornnccf, 

J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 
Calzetti, D. ct al. 2005, ApJ, 633, 871 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Chariot, S. ., & Fall, S. M. 2000, ApJ, 539, 718 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 
Code, A. D., & Welch, G. A. 1982, ApJ, 256, 1 
Condon, J. J. 1992, ARA&A, 30, 575 

Consclice, C. J., Gallagher, J. S., Calzctti, D., Homeier, N., & 

Kinney, A. 2000, AJ, 119, 79 
Dale, D. A. et al. 2005, ApJ, 633, 857 
Dale, D. A., & Helou, G. 2002, ApJ, 576, 159 

de Vaucouleurs, G., de Vaucoulcurs, A., Corwin, H. G., Buta, R. J., 
Paturel, G., & Fouque, P. 1992, VizieR Online Data Catalog, 
7137, 

Devereux, N. A., Jacoby, G., & Ciardullo, R. 1995, AJ, 110, 1115 
Devriendt, J. E. G., Guiderdoni, B., & Sadat, R. 1999, A&A, 350, 
381 

Draine ct al. 2006, ApJ (in preparation) 

Dunne, L., Ealcs, S., Edmunds, M., Ivison, R., Alexander, P., & 

Clements, D. L. 2000, MNRAS, 315, 115 
Dunne, L., & Eales, S. A. 2001, MNRAS, 327, 697 
Elvis, M., & van Speybroeck, L. 1982, ApJ, 257, L51 
Engelbracht, C. W. et al. 2004, ApJS, 154, 248 
Engelbracht, C. W., Gordon, K. D., Rieke, G. H., Werner, M. W., 

Dale, D. A., & Latter, W. B. 2005, ApJ, 628, L29 
Fabbiano, G. 1988, ApJ, 325, 544 
Fazio, G. G. et al. 2004, ApJS, 154, 10 
Ferland, G. J. 1980, PASP, 92, 596 

Filippenko, A. V., & Sargent, W. L. W. 1988, ApJ, 324, 134 
Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Flores, H. et al. 1999, ApJ, 517, 148 

Forster Schreiber, N. M., Roussel, H., Sauvage, M., & 

Charmandaris, V. 2004, A&A, 419, 501 
Franccschini, A., Aussel, H., Cesarsky, C. J., Elbaz, D., & Fadda, 

D. 2001, A&A, 378, 1 
Freedman, W. L. et al. 1994, ApJ, 427, 628 
Garnett, D. R., & Shields, G. A. 1987, ApJ, 317, 82 
Genzel, R., & Cesarsky, C. J. 2000, ARA&A, 38, 761 
Gil de Paz, A., Aragon-Salamanca, A., Gallego, J., Alonso-Herrero, 

A., Zamorano, J., & Kauffmann, G. 2000, MNRAS, 316, 357 
Gil dc Paz, A., & Madore, B. F. 2002, AJ, 123, 1864 
Goldadcr, J. D., Mcurcr, G., Heckman, T. M., Seibcrt, M., Sanders, 

D. B., Calzetti, D., & Steidel, C. C. 2002, ApJ, 568, 651 



Gordon, K. D., Clayton, G. C, Witt, A. N., & Missclt, K. A. 2000, 

ApJ, 533, 236 
Gordon, K. D. et al. 2004, ApJS, 154, 215 
— . 2005, PASP, 117, 503 
Gordon ct al. 2006, ApJ (submitted) 
Gottesman, S. T., & Weliachew, L. 1975, ApJ, 195, 23 
Greenawalt, B., Waltcrbos, R. A. M., Thilker, D., & Hoopes, C. G. 

1998, ApJ, 506, 135 
Heavens, A., Panter, B., Jimenez, R., & Dunlop, J. 2004, Nature, 

428, 625 

Heckman, T. M. 1980, A&A, 87, 152 
Helou, G. et al. 2004, ApJS, 154, 253 
Hill, J. K. et al. 1992, ApJ, 395, L37 
— . 1995, ApJ, 438, 181 
Hinz, J. L. et al. 2004, ApJS, 154, 259 

Hippclein, H, Haas, M., Tuffs, R. J., Lemke, D., Stickel, M., Klaas, 

U., & Volk, H. J. 2003, A&A, 407, 137 
Ho, L. C, Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 

112, 315 

H0g, E. et al. 2000, A&A, 357, 367 

Holwcrda, B. W., Gonzalez, R. A., van dcr Kruit, P. C, & Allen, 

R. J. 2005, A&A, 444, 109 
Hopkins, A. M. 2004, ApJ, 615, 209 

Hopkins, A. M., Connolly, A. J., Haarsma, D. B., & Cram, L. E. 

2001, AJ, 122, 288 
Houck, J. R. et al. 2004, ApJS, 154, 211 

Jarrett, T. H., Chester, T., Cutri, R., Schneider, S., Skrutskie, M., 

& Huchra, J. P. 2000, AJ, 119, 2498 
Karachentsev, I. D. et al. 2002, A&A, 383, 125 
Kauffmann, G. et al. 2003, MNRAS, 341, 33 

Kaufman, M., Bash, F. N., Kennicutt, R. C, & Hodge, P. W. 1987, 

ApJ, 319, 61 
Kennicutt, R. C. 1998, ARA&A, 36, 189 
Kennicutt, R. C. ct al. 2003, PASP, 115, 928 
Kennicutt ct al. 2006, ApJ (in press) 

Kewley, L. J., Geller, M. J., Jansen, R. A., & Dopita, M. A. 2002, 

AJ, 124, 3135 
Kleinmann, D. E., & Low, F. J. 1970, ApJ, 161, L203 
Kong, X., Chariot, S., Brinchmann, J., & Fall, S. M. 2004, MNRAS, 

349, 769 

Kong, X. et al. 2000, AJ, 119, 2745 

La Parola, V., Fabbiano, G., Elvis, M., Nicastro, F., Kim, D. W., 

& Peres, G. 2004, ApJ, 601, 831 
Laor, A., & Draine, B. T. 1993, ApJ, 402, 441 
Leger, A., & Puget, J. L. 1984, A&A, 137, L5 
Li, A., & Draine, B. T. 2001, ApJ, 554, 778 
Lin, W. et al. 2003, AJ, 126, 1286 
Martin, D. C. et al. 2005, ApJ, 619, LI 

Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64 
Morrissey, P. et al. 2005, ApJ, 619, L7 

Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and 
active galactic nuclei (Research supported by the University 
of California, John Simon Guggenheim Memorial Foundation, 
University of Minnesota, ct al. Mill Valley, CA, University 
Science Books, 1989, 422 p.) 

Perez-Gonzalez, P. G., Gallego, J., Zamorano, J., Alonso-Hcrrcro, 
A., de Paz, A. G., & Aragon-Salamanca, A. 2003a, ApJ, 587, 
L27 

Perez-Gonzalez, P. G., Gil de Paz, A., Zamorano, J., Gallego, J., 
Alonso-Herrero, A., & Aragon-Salamanca, A. 2003b, MNRAS, 
338, 508 

— . 2003c, MNRAS, 338, 525 

Perez-Gonzalez, P. G. et al. 2005, ApJ, 630, 82 

Perez-Gonzalez, P. G., Zamorano, J., Gallego, J., Aragon- 
Salamanca, A., & Gil de Paz, A. 2003d, ApJ, 591, 827 



22 



Perez-Gonzalez et al. 



Page, M. J., Breeveld, A. A., Soria, R., Wu, K., Branduardi- 
Raymont, G., Mason, K. O., Starling, R. L. C, & Zane, S. 2003, 
A&A, 400, 145 

Pahre, M. A., Ashby, M. L. N., Fazio, G. G., & Willner, S. P. 2004, 
ApJS, 154, 235 

Papovich, C, Dickinson, M., & Ferguson, H. C. 2001, ApJ, 559, 
620 

Peeters, E., Spoon, H. W. W., & Tielens, A. G. G. M. 2004, ApJ, 
613, 986 

Peimbert, M., & Torres-Peimbert, S. 1981, ApJ, 245, 845 

Petre, R., Mushotzky, R. F., Serlemitsos, P. J., Jahoda, K., & 

Marshall, F. E. 1993, ApJ, 418, 644 
Popescu, C. C, & Tuffs, R. J. 2003, A&A, 410, L21 
Popescu, C C et al. 2005, ApJ, 619, L75 

Popescu, C. C, Tuffs, R. J., Volk, H. J., Pierini, D., & Madorc, 

B. F. 2002, ApJ, 567, 221 
Prescott et al. 2006, ApJ (in press) 
Puget, J. L., & Leger, A. 1989, ARA&A, 27, 161 
Reach, W. T. et al. 2005, PASP, 117, 978 
Regan, M. W. et al. 2004, ApJS, 154, 204 

Reichen, M., Kaufman, M., Blecha, A., Golay, M., & Huguenin, D. 

1994, A&AS, 106, 523 
Rice, W., Lonsdale, C. J., Soifer, B. T., Neugebauer, G., Koplan, 

E. L., Lloyd, L. A., dc Jong, T., & Habing, H. J. 1988, ApJS, 

68, 91 

Rieke, G. H., & Low, F. J. 1972, ApJ, 176, L95 

Rieke, G. H. et al. 2004, ApJS, 154, 25 

Rots, A. H., & Shane, W. W. 1975, A&A, 45, 25 

Roussel, H., Gil de Paz, A., Seibert, M., Helou, G., Madorc, B. F., 

& Martin, C. 2005, ApJ, 632, 227 
Rowan-Robinson, M. 2001, ApJ, 549, 745 
Sage, L. J., & Westpfahl, D. J. 1991, A&A, 242, 371 
Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 



Saunders, W., Rowan- Robinson, M., Lawrence, A., Efstathiou, G., 
Kaiser, N., Ellis, R. S., & Frenk, C. S. 1990, MNRAS, 242, 318 
Sauvage, M., & Thuan, T. X. 1992, ApJ, 396, L69 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Shuder, J. M., & Osterbrock, D. E. 1981, ApJ, 250, 55 
Smith, J. D. T. et al. 2004, ApJS, 154, 199 
Solinger, A., Morrison, P., & Markert, T. 1977, ApJ, 211, 707 
Spoon, H. W. W. et al. 2004, ApJS, 154, 184 
Stauffcr, J. R., & Bothun, G. D. 1984, AJ, 89, 1702 
Stickel, M., Barnes, D., & Krause, O. 2005, A&A, 443, 373 
Sturm, E., Lutz, D., Tran, D., Feuchtgruber, H., Genzel, R., Kunze, 
D., Moorwood, A. F. M., & Thornlcy, M. D. 2000, A&A, 358, 
481 

Sullivan, M., Mobasher, B., Chan, B., Cram, L., Ellis, R., Treyer, 

M., & Hopkins, A. 2001, ApJ, 558, 72 
Walterbos, R. A. M., & Greenawalt, B. 1996, ApJ, 460, 696 
Weingartner, J. C, & Draine, B. T. 2001, ApJ, 548, 296 
Werner, M. W. et al. 2004a, ApJS, 154, 1 

Werner, M. W., Uchida, K. I., Sellgren, K., Marengo, M., Gordon, 
K. D., Morris, P. W., Houck, J. R., & Stansberry, J. A. 2004b, 
ApJS, 154, 309 

Williams, G. G., Olszewski, E., Lesser, M. P., & Burge, J. H. 2004, 
in Ground-based Instrumentation for Astronomy. Edited by Alan 
F. M. Moorwood and lye Masanori. Proceedings of the SPIE, 
Volume 5492, pp. 787-798 (2004)., 787-798 

Witt, A. N., & Gordon, K. D. 2000, ApJ, 528, 799 

Worthey, G. 1994, ApJS, 95, 107 

Wu, C.-C, Gallagher, J. S., Peck, M., Fabcr, S. M., & Tinsley, 

B. M. 1980, ApJ, 237, 290 
Youngblood, A. J., & Hunter, D. A. 1999, ApJ, 519, 55 
Zaritsky, D., Kennicutt, R. C, & Huchra, J. P. 1994, ApJ, 420, 87 



